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This  dissertation  presents  methodology  for  physical  charge- 
based  modeling  of  a lateral  insulated-gate-bipolar  transistor  (LIGBT) 
for  integrated  circuit  computer  design.  The  steady-state  and  transient 
carrier  transport  problems  in  the  LIGBT  are  formulated  and  solved. 
The  effects  of  the  two-dimensional  current  analysis  in  the  lateral  PNP 
BJT/LIGBT  subcell  is  characterized,  and  the  current-induced-space- 
charge  region  in  the  LDMOST/LIGBT  subcell  is  modeled.  These 
effects,  which  are  not  presented  in  conventional  equivalent- 
(sub)circuit  models,  are  physically  and  sometimes  semi-numerically 
accounted  for  in  our  model.  Improvements  on  both  subcell  models  of 
the  LIGBT  have  resulted  in  accurate  simulation  results  without 
parameter  optimization.  The  Kirk  effect  is  modeled  and  added  in  our 
LIGBT  model.  Two-dimensional  numerical  device  simulations  were 
used  extensively  to  study  the  effects  and  to  aid  the  model 


development.  The  LIGBT  model  is  implemented  in  the  Saber  circuit 
simulator,  by  which  the  model  equations  are  solved  semi-numerically 
within  the  nodal  analysis  framework  of  the  Saber.  The  parameter 
extraction  method  for  Saber  simulation  is  developed.  Our  new  model  is 
verified  by  measurements  of  the  test  devices,  using  simulations  with 
model  parameters  extracted  from  static  and  dynamic  measurements 
without  any  optimization  to  fit  the  measured  data.  The  LIGBT  model 
in  Saber  provides  the  capability  for  mixed-mode  device/circuit 
simulation  and  hence  can  facilitate  computer-aided  optimal  device/ 
circuit  design  of  IGBTs. 


VI 


CHAPTER  1 
INTRODUCTION 


Recently,  advanced  process  technologies  of  high-voltage  (HV) 
devices  have  enabled  the  integration  of  HV  devices  with  low-voltage 
control  circuitry  on  the  same  chip.  These  technologies  usually  mix  con- 
ventional bipolar  and  MOS  processes  and  utilize  advanced  VLSI 
technology  to  create  a very  promising  area  in  the  IC  industry.  The 
application  areas  include  high-voltage  displays,  motor  controllers,  reg- 
ulators, TV  circuits,  telecommunication,  switches,  and  interface-chips 
[BEC85,  RUM85,  ADL84]. 

The  lateral  insulated-gate  bipolar  transistor  (LIGBT)  is  one 
of  the  MOS-controlled  HV  devices  for  power  IC  applications.  The  MOS 
channel  in  the  LIGBT  is  usually  formed  by  double-diffusion  process; 
sequential  diffusion  of  p-  and  n-type  impurities  in  the  epi  region  yields 
the  MOS  channel.  The  double-diffusion  process  is  widely  used  for  HV 
devices  because  of  easier  process  steps  to  achieve  short  channels  while 
realizing  high  breakdown  voltages. 

The  LIGBT  combines  the  advantages  of  current-carrying 
capability  and  the  conductivity  modulation  of  the  IGBT  with  high 
input  impedance  and  the  high  breakdown  voltage  of  RESURF 
LDMOST.  The  on-resistance  is  about  an  order  of  magnitude  lower  than 
the  LDMOST  transistor  of  equivalent  dimensions  because  the  base 
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region  is  conductivity  modulated.  The  cross-sectional  view  of  the 
LIGHT  device  is  shown  in  Fig.  1.1.  This  device  was  fabricated  with 
dielectrically  isolated  bipolar-CMOS-DMOS  (BCDMOS)  integrated-cir- 
cuit technology  developed  at  AT&T  Bell  Laboratories  for  HV 
applications  [LU88].  Dielectric  isolation  (DI)  effectively  minimize  par- 
asitic cross  coupling  between  circuit  elements  and  to  reduce  the 
separation  space  between  HV  components  compared  to  that  which  can 
be  achieved  with  conventional  junction  isolation  techniques. 

Along  with  the  matured  process  technology  of  the  LIGHT  and 
integrated  circuits,  the  accurate  CAD  simulation  tool  for  the  LIGBT 
should  be  developed  to  support  optimal  device  and  circuit  design.  The 
circuit  simulation  for  the  LIGBT  still  relies  on  SPICE-based  device 
models  due  to  the  unavailability  of  accurate  the  LIGBT  device  model 
[HEF93].  The  SPICE-based  models  are  unable  to  predict  large  over- 
shoot voltages  and  currents  and  other  switching  effects  critical  in  the 
design  of  power  converters.  Therefore,  the  development  of  physics- 
based  model  for  the  LIGBT  and  implementation  into  the  circuit  simu- 
lation tool,  e.g.,  SABER,  are  needed.  The  advantages  of  a physics- 
based  model  are  the  accuracy  of  simulation  results  and  the  fundamen- 
tal parameters  which  are  available  from  layout  and  process 
information. 

The  purpose  of  this  dissertation  is  to  develop  and  implement 
a physical  model  for  the  LIGBT.  The  effects  of  the  current-induced- 
space-charge  region  for  the  LDMOST/LIGBT  subcell  and  two-dimen- 
sional current  flow  for  a lateral  PNP  BJT/LIGBT  subcell  are  modeled. 
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Fig.  1.1  The  cross  sectional  view  of  unit-cell  of  n-channel  LIGBT. 

The  series  resistor  load  and  gate  drive  test  circuitry  corre- 
spond to  a transient  turn-off  from  steady-on-state. 
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and  the  Kirk  effect  is  characterized.  The  model  parameter  extraction 
method  is  developed  and  the  physics-based  LIGBT  model  is  imple- 
mented in  the  Saber  circuit  simulator.  The  major  contributions  made 
in  this  dissertation  are 

(1)  the  development  of  the  generalized  methodology  for  phys- 
ics-based modeling  of  IGBTs  based  on  two-dimensional 
numerical  device  simulations,  and  the  methodology  can  be 
useful  to  model  other  power  devices; 

(2)  the  formulation  and  solution  of  the  steady-state  and  tran- 
sient carrier  transport  problems  in  the  LIGBT,  and 
improvement  on  both  LDMOST/LIGBT  subcell  and  lateral 
PNP  BJT/LIGBT  subcell  models; 

(3)  the  development  of  a parameter  extraction  method  for  the 
LIGBT; 

(4)  the  implementation  of  the  physics-based  LIGBT  model 
into  the  Saber  circuit  simulator  and  the  model  verification 
for  various  external  circuit  conditions. 

In  Chapter  2,  the  generalized  methodology  for  physics-based 
modeling  of  LIGBT  is  presented.  The  generalized  methodology  implies 
not  only  identification  of  each  basic  semiconductor  region,  but  also 
application  of  semiconductor  equations  which  represent  each  semicon- 
ductor region.  The  study  gives  physical  insight  based  on  extensive 
two-dimensional  numerical  device  simulations  and  describes  proper 
physics-based  modeling  techniques  of  IGBT.  The  generalized  method- 
ology for  physics-based  modeling  developed  in  Chapter  2 can  be  useful 
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to  model  the  power  devices  with  wide  and  lightly  doped  regions  such  as 
the  power  bipolar  transistor  and  the  PIN  diode. 

In  Chapter  3,  a physical  model  for  the  LIGBT  is  described. 
Improvements  on  both  LDMOST/LIGBT  and  lateral  PNP  BJT/LIGBT 
subcell  models  have  resulted  in  accurate  simulation  results  without 
parameter  optimization.  In  the  LDMOST/LIGBT  subcell,  we  have  mod- 
eled the  current-induced-space-charge  region  shown  to  be  of  great 
importance  for  improving  the  accuracy  of  the  LDMOST  model.  In  the 
lateral  PNP  BJT/LIGBT  subcell,  we  have  derived  analytical  geometri- 
cal factors  based  on  the  ambipolar  transport  equation.  The  lateral  and 
quasi-vertical  geometrical  factors  enhance  the  analysis  of  the  two- 
dimensional  anode  current,  and  thus,  result  in  accurate  modeling  for 
steady-state  and  transient  conditions  for  the  LIGBT.  The  implementa- 
tion of  the  Kirk-effect  into  the  LIGBT  model  also  improved  the  model 
accuracy  of  steady-state  and  transient  conditions.  PISCES  simulations 
and  measurements  of  test  devices  support  the  LIGBT  model,  which  is 
demonstrated  in  steady-state  and  turn-off  transient  Saber 
simulations. 

In  Chapter  4,  the  parameter  extraction  method  for  the  LIGBT 
from  terminal  electrical  measurements  has  been  developed  and 
described.  Unlike  microelectronic  devices,  the  dynamic  characteristic 
must  be  used  to  extract  model  parameters  in  the  LIGBT.  This  is 
because  the  LIGBT  has  the  characteristics  of  the  internal  LDMOST 
and  lateral  PNP  bipolar  transistor  subcells.  The  dynamic  waveforms 
contain  many  features  that  isolate  different  model  parameters.  The 
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model  parameters  are  obtained  sequentially  by  selecting  features  of 
the  device  characteristics  that  isolate  parameters,  and  using  the 
parameters  obtained  from  previous  extraction  steps  to  calculate  the 
model  parameters  from  the  next  measured  characteristics. 

In  Chapter  5,  the  methodology  for  implementing  the  physical 
model  for  the  LIGBT  into  the  Saber  circuit  simulator  is  described.  The 
Saber  is  a very  powerful  modeling  tool  for  power  devices  which  have 
widely  varying  cell  structures  even  within  one  category  of  device.  The 
numerical  machine  in  the  Saber  is  transparent  to  the  user,  and  one 
template  containing  the  model  equations  is  the  only  implementation 
needed.  The  techniques  for  implementing  the  LIGBT  model  equations 
into  the  Saber  circuit  simulator  are  essential  for  accurate  simulations. 
Furthermore,  the  techniques  demonstrated  for  modeling  the  LIGBT 
with  the  Saber  circuit  simulator  are  generally  applicable  to  modeling 
other  power  semiconductor  devices.  The  LIGBT  model  is  verified  for 
various  external  circuit  conditions  and  for  the  full  range  of  steady-sate 
and  dynamic  conditions  in  which  the  LIGBT  is  intended  to  be  oper- 
ated. The  LIGBT  Saber  model  also  is  shown  to  be  suitable  for 
simulating  the  behavior  of  the  LIGBT  for  general  purpose  external  cir- 
cuit conditions  and  designing  the  various  power  circuits. 

In  Chapter  6,  the  main  accomplishments  of  this  dissertation 
are  summarized,  and  suggestions  for  further  research  are  discussed. 

Appendix,  related  to  Chapter  5,  lists  the  source  code  of  the 


LIGBT  Saber  model  template. 


CHAPTER  2 

GENERALIZED  METHODOLOGY  FOR  PHYSICS-BASED 
MODELING  OF  IGBTS 

2.1  Introduction 

In  the  1980s,  the  Insulated  Gate  Bipolar  Transistor  (IGBT) 
was  introduced  to  overcome  the  high  on-state  power  loss  of  the  power 
MOSFET  while  maintaining  the  simple  gate  drive  requirements  of  the 
device.  Because  of  the  special  structure  of  the  IGBTs,  the  standard 
models  designed  for  microelectronic  devices  are  not  suitable  to 
describe  the  dynamic  behavior  of  the  IGBTs.  Unlike  low  voltage  IC 
devices,  IGBTs  have  wide  and  lightly  doped  base  regions  to  support 
high-voltage  operation.  This  lightly  doped  high-level  injection  region 
significantly  impacts  the  electrical  behavior  of  the  IGBTs,  and  is  diffi- 
cult to  model  because  its  distributed  phenomenon  must  be  described 
by  the  ambipolar  transport  equation.  This  is  the  main  reason  that 
computer-aided-design  for  IGBTs  is  hardly  used  in  the  power  elec- 
tronic field  today.  For  the  design  and  analysis  of  power  ASICS,  we 
need  to  develop  accurate  models  for  the  IGBT  and  implement  the  mod- 
els into  the  circuit  simulators  like  SABER  and  SPICE. 

In  this  chapter,  we  present  the  generalized  methodology  of 
physics-based  modeling  for  IGBTs.  It  starts  from  identification  of  the 
basic  semiconductor  regions  in  IGBT’s  mode  of  operation.  The  basic 
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semiconductor  regions  can  be  categorized  into  five  parts  in  IGBTs:  the 
anode-base  and  the  base-cathode  p/n  junction  region,  the  anode  and 
cathode  terminal  doping  region,  the  high-level  injection  base  ambipo- 
lar  region,  the  current-induced-space-charge  region  at  the  edge  of  the 
LDMOST  channel,  and  the  channel  region.  Each  basic  semiconductor 
region  can  be  identified  based  on  the  two-dimensional  PISCES  simula- 
tion. As  a matter  of  fact,  numerous  works  about  regional  modeling  of 
IGBTs  [GOE94,  ALL94,  MA94]  have  been  proposed  while  no  works 
focused  on  the  generalized  methodology  of  modeling  for  IGBTs.  The 
generalized  methodology  implies  not  only  identification  of  each  basic 
semiconductor  region  but  also  application  of  semiconductor  equations 
in  the  basic  semiconductor  regions.  We  have  four  major  semiconductor 
equations:  the  current  density  equations,  the  current  continuity  equa- 
tions, Poisson’s  equation,  and  Boltzmann’s  relation.  Different  physical 
phenomena  can  be  observed  in  different  regions  of  the  semiconductor. 
Therefore,  some  semiconductor  equations  can  be  neglected  in  the  mod- 
eling of  each  semiconductor  region.  This  can  make  the  modeling  of 
power  devices  simple  and  easy  enough  for  circuit  design.  This  general- 
ized methodology  can  be  applied  to  model  the  power  devices  that  have 
wide  and  lightly  doped  base  region  such  as  the  power  bipolar  transis- 
tor and  the  PIN  diode. 

Figure  2.1  summarizes  the  unified  view  of  the  generalized 
methodology  for  physics-based  modeling  for  the  IGBTs.  The  first  step 
is  to  collect  the  device  information,  i.e.,  its  geometry,  layout,  and 
dopant  profile.  The  layout  and  processing  information  input  into  the 
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Fig. 2.1  Flowchart  of  the  generalized  methodology  of  physics- 

based  modeling  of  IGBTs. 
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PISCES  simulator  to  provide  insight  into  the  device  operation.  Based 
on  extensive  PISCES  simulations,  we  can  identify  the  basic  semicon- 
ductor regions  in  the  desired  modes  of  operation.  Besides  identifying 
the  important  semiconductor  regions  in  the  device,  the  two-dimen- 
sional PISCES  simulation  can  be  useful  in  determining  the  sizes  and 
the  shapes  of  each  semiconductor  region.  Then,  we  can  apply  the  semi- 
conductor equations  which  represent  each  semiconductor  region.  In 
some  regions,  only  one  main  semiconductor  equation  can  be  applied 
because  one  main  physical  phenomenon  is  observed.  By  solving  the 
equations  in  each  semiconductor  region  and  linking  them  with  the 
boundary  conditions,  we  can  make  a regional  model  for  the  IGBT. 


2.2  Identification  of  Basic  Semiconductor  Regions 

The  identification  of  the  basic  semiconductor  regions  is  an 
important  part  of  the  physics-based  modeling  of  power  devices.  As 
stated  in  Section  2.1,  we  have  five  different  basic  semiconductor 
regions  in  IGBTs.  In  this  section,  the  identification  of  the  basic  semi- 
conductor regions  during  the  IGBT’s  operation  is  described.  This  study 
is  based  on  extensive  simulations  of  the  IGBT  using  the  two-dimen- 
sional device  simulator  PISCES,  which  utilizes  finite-element 
approximations  and  a Gummel-Newton  numerical  iteration  method  to 
solve  the  discretized  semiconductor  equations  defining  the  carrier 
transport  [COU88].  Doping  density  profiles  and  structural  dimensions 
for  our  test  devices  were  obtained  from  AT&T  Bell  laboratories.  The 
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simulations  reveal  the  internal  mechanisms  of  IGBT  such  as  different 
collecting  regions  of  the  holes  and  the  electrons,  and  they  show  that 
the  current-induced-space-charge  region  forms  at  the  edge  of  the 
LDMOST  channel  of  the  LIGBT  because  of  its  similar  structure  to  the 
RESURF  LDMOST.  Only  a depletion  region  is  observed  in  the  base 
area  of  the  VIGBT  [HEF88].  Hence,  the  current  transport  for  the 
LIGBT  and  the  VIGBT  are  very  different  from  each  other. 

The  PISCES  input  file  for  the  LIGBT  is  listed  in  the  Appendix 
A.  The  geometrical  dimensions  and  the  doping  density  profiles  for  each 
region  are  specified  in  the  input  file.  This  device  was  fabricated  with  a 
dielectric  isolation  (DI)  Bipolar-CMOS-DMOS  technology  developed  at 
AT&T  Bell  Laboratories  [KOS93].  Figure  2.2  shows  the  cross-section  of 
the  LIGBT  generated  by  the  PISCES  input  file.  From  the  doping  pro- 
files, the  LIGBT  has  five  different  doped  regions:  anode,  deep  p'*’,  p’ 
body,  source,  and  the  n'  uniform  doped  base,  and  it  has  three  elec- 
trodes: anode,  cathode,  and  gate. 

As  shown  in  Fig. 2. 2,  the  LIGBT  has  two  subcells  indicating 
an  LDMOS  gate  structure  with  the  bipolar  current  conduction  of  the 
lateral  PNP  BJT.  In  this  structure,  the  current  is  negligible  when  a 
negative  voltage  is  applied  to  the  anode  with  respect  to  the  cathode 
because  the  anode-base  junction  will  become  reverse-biased.  This  pro- 
vides the  LIGBT  with  its  reverse  blocking  capability.  When  a positive 
voltage  is  applied  to  the  anode  terminal  with  the  gate  short-circuited 
to  the  cathode  terminal,  the  base-cathode  junction  becomes  reverse- 
biased,  and  the  device  operates  in  its  forward  blocking  mode.  At  posi- 
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CROSS-SECTION  OF  UGBT 


CATHODE  GATE  ANODE 


Fig. 2. 2 Cross-section  of  the  AT&T  LIGBT  (not  scaled).  Five  differ- 
ent doped  regions  and  three  electrodes:  anode,  cathode, 
and  gate. 
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tive  anode  terminal  voltages,  if  a positive  gate  bias  is  applied  of 
sufficient  magnitude  to  invert  the  surface  of  the  p'  body  region  under 
the  gate  (channel),  the  LIGBT  operates  in  its  forward  conducting  state 
because  electrons  can  now  flow  from  the  n'*’  source  region  to  the  n'  epi 
base  region.  This  is  well  shown  in  Fig.  2.3(a).  The  electron  current  vec- 
tor appears  only  in  the  source  and  the  channel  regions.  In  this  forward 
conducting  state,  the  anode-base  junction  is  forward-biased  and  the 
anode  p"*"  region  injects  holes  into  the  n*  epi  base  region  shown  in  Fig. 
2.3(b).  The  hole  current  vector  predominates  throughout  the  entire 
LIGBT  structure  (except  in  the  source  region)  in  the  forward  conduct- 
ing mode. 

The  LIGBT  has  a LDMOST/LIGBT  subcell  and  a lateral  PNP 
BJT/LIGBT  subcell.  When  the  forward  bias  (anode  voltage)  is 
increased,  the  injected  hole  concentration  increases  until  it  exceeds 
the  background  doping  level  of  the  n’  base  region.  When  the  excess 
carrier  density  exceeds  the  background  doping  level  by  several  orders 
of  magnitude,  as  is  the  case  in  the  injected  part  of  the  base  region,  the 
high-injection  condition  is  satisfied.  The  transport  of  electrons  and  the 
transport  of  holes  are  coupled  by  the  electric  field  in  the  drift  region, 
and  they  cannot  be  treated  separately.  Consequently,  the  basic  equa- 
tion governing  the  carrier  transport  is  the  ambipolar  transport 
equation.  Figure  2.4(b)  shows  the  one-dimensional  carrier  concentra- 
tion at  y=2pm,  and  it  clearly  shows  that  the  doping  concentration  of 
the  injected  carriers  are  much  higher  than  the  base  background  doping 
concentration  (=2xl0^‘^cm'^)  and  the  electron  and  hole  concentrations 


Distance  (Microns)  Distance  (Microns) 
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(a) 


LIGBT  Electron  Current  Vector  at  Va=10  V Vg=10  V 
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LIGBT  Hole  Current  Vector  at  Va=10  V Vg=10  V 
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Fig. 2. 3 (a)  Electron  and  (b)  Hole  current  vectors  of  LIGBT  at  Va=10 
V and  Vg=10  V 
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(a) 

LIGBT  Depletion  region  Va=10  V Vg=10  V 


(b) 

LIGBT  one-dimensional  carrier  concentration  at  y=6pm. 


Fig.2.4  (a)  Depletion  region  in  the  LIGBT  at  Va=10  V and  Vg=10  V. 

(b)  One  dimensional  carrier  distribution  in  the  base  at  y=6M,m. 
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in  the  base  region  is  equal  in  = p).  From  this,  we  know  the  base  region 
is  a high-level-injected  ambipolar  transport  region. 

The  LDMOST/LIGBT  subcell  has  the  structure  of  the 
RESURF  LDMOST.  The  RESURF  LDMOST  has  the  n'epi  layer  which 
functions  as  a drift  region.  The  thickness  and  doping  density  of  the  n' 
epi  layer  are  controlled  to  produce  a low  surface  electric  field  to  pro- 
vide high  breakdown  voltage  [APP79].  At  high  currents,  excess 
electrons  accumulate  in  the  drift  region  due  to  drift-velocity  satura- 
tion. The  excess  electrons  produce  space-charge-limited  current  flow 
and  compensate  for  the  widening  of  the  depletion  region  of  base-collec- 
tor junction  in  the  lateral  PNP  BJT/LIGBT  subcell.  Figure  2.5(a) 
shows  the  equipotential  distribution  of  two-dimensional  PISCES  simu- 
lation of  the  LIGBT  mode  of  operation  at  = 10  V and  = 10  V. 
The  dotted  line  indicates  the  edge  of  the  base-collector  depletion 
region.  As  shown  by  PISCES  simulation,  the  equipotential  contour  not 
only  is  distributed  along  the  base-collector  depletion  region  but  also  is 
distributed  along  the  quasi-neutral  base  region  near  the  channel.  The 
equipotential  contour  along  the  quasi-neutral  base  has  a cone-shaped, 
cross-sectional  area.  The  cone-shaped  area  is  a current-induced-space- 
charge  (CISC)  region  which  is  incorporated  into  the  base  current.  Fig- 
ure 2.5(b)  shows  the  PISCES  simulation  of  a one-dimensional  electric 
field  distribution  at  the  surface  (y=0.001|im).  In  Fig.2.5(a),  the  physi- 
cal location  of  the  drain  of  LDMOST/LIGBT  subcell  is  shown. 

The  simulated  PISCES  current  flowlines  within  the  LIGBT 
device  is  shown  in  Fig. 2. 6.  This  simulated  current  flowlines  reflect 
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(a) 

LIGHT  Potential  Distribution  Va=10  V and  Vg=10  V 


Distance  (Microns) 


(b) 

LIGHT  one-dimensional  Electric  Field  Distribution  at  y=0|i,m 


Fig. 2. 5 (a)  Equipotential  contour  from  two-dimensional  PISCES 
simulation  of  the  LIGHT  mode  of  operation. 

(b)  One-dimensional  electric  field  distribution  at  the 
surface. 
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(a) 


LIGBT  current  flow  at  Va=10  V and  Vg=10  V 


Fig. 2. 6 Anode  current  flowlines  from  two-dimensional  PISCES 
simulation  at  Vgs=10  V and  Vak  = 10  V.  The  flowlines  are 
plotted  such  that  the  equal  currents  (A/ = 0.2  x 10~^A/pm) 
flow  between  any  two  adjacent  flowlines  (a)  LIGBT 
(scaled)  (b)  SOI  LIGBT  (not  scaled). 
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operation  below  the  onset  of  latchup.  Electrons  injected  (predomi- 
nantly laterally)  from  the  source  drive  the  lateral  PNP  BJT/LIGBT 
base  harder,  producing  additional  forward  bias  on  the  anode  junction 
and  additional  hole  injection  into  the  n*  base.  This  enhances  the  con- 
ductivity modulation  of  the  lateral  PNP  BJT/LIGBT  base  region, 
including  the  adjacent  depletion  region  at  the  lateral  PNP  BJT/LIGBT 
base-collector  junction.  Normal  LIGBT  operation  is  evident;  the 
LDMOST  current  through  the  channel  supplies  electrons  to  the  base  of 
the  lateral  PNP/BJT  subcell  which  then  supports,  through  recombina- 
tion, the  hole  transport  current.  The  LIGBT  anode  current  is  the  sum 
of  the  LDMOST  channel  current  and  the  collector  current  of  the  lat- 
eral PNP  BJT/LIGBT  subcell.  Note  that  all  the  current  flowline  are 
not  lateral,  but  some  have  quasi-vertical  flowlines.  The  SOI  LIGBT 
[DIS93],  in  Fig.  2.6(b),  has  predominantly  lateral  flowlines.  This  is 
mainly  because  the  SOI  LIGBT  has  only  2pm  of  n"  epi  base  region.  The 
geometrical  parameters  such  as  base  width,  junction  depth,  and  base 
length  influence  the  current  flowlines.  Hence,  we  need  to  model  two- 
dimensional  analysis  of  current  flowlines  which  are  dependent  on  the 
geometrical  input  parameters. 

Now  that  the  identiflcation  of  the  basic  semiconductor  regions 
in  the  LIGBT  is  completed,  and  the  circuit  models  for  each  building 
blocks  can  be  defined.  Figure  2.7  shows  the  block  diagram  of  the  basic 
semiconductor  regions  in  the  LIGBT  mode  of  operation.  It  consists  of 
seven  blocks,  and  each  block  represents  a basic  semiconductor  region. 
Among  the  basic  semiconductor  regions,  the  high-level  injection  ambi- 
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Fig. 2. 7 Block  diagram  of  basic  semiconductor  regions  for  the 
LIGBT  mode  of  operation 
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polar  base  region  is  the  most  important  region  because  it  determines 
the  electrical  behavior  of  the  LIGBT.  In  next  section,  the  application  of 
semiconductor  equations  to  the  basic  semiconductor  region  is 
illustrated. 

The  PISCES  input  file,  which  details  the  geometry  input  and 
processing  information,  is  listed  in  Appendix  A. 

2.3  Application  of  Semiconductor  Equations 

As  stated  in  Section  2.2,  the  LIGBT  exhibits  distinctive  basic 
semiconductor  regions,  such  as  the  lightly  doped  cm’^)  base 

region  for  blocking  high  voltages  where  most  of  the  excess  carriers  are 
stored,  the  more  heavily  doped  (-10^®  cm'^)  end  regions  for  the  termi- 
nals, and  the  lightly  doped  (-10^®  cm‘^)  semiconductor  surface  in  an 
MOS  structure.  The  semiconductor  equations  - the  current  density 
equations,  the  current  continuity  equations,  Poisson’s  equation,  and 
Boltzmann’s  equation-  do  not  have  the  same  importance  in  the  differ- 
ent regions  of  the  IGBTs.  Some  of  the  semiconductor  equations  can  be 
neglected  in  modeling  each  semiconductor  region.  In  Table  2.1,  the 
basic  semiconductor  regions  with  identification  criteria  and  the  semi- 
conductor equations  which  fit  each  region  are  listed.  Except  the 
channel  region,  all  other  regions  are  related  to  the  behavior  of  the 
bipolar  power  semiconductor  devices. 
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Table  2.1  The  main  types  of  behavior  inside  IGBTs 


Region 

Identification  criteria 

Equations 

P/N  junction 
region 

Seperation  of  quasi-fermi  level  in  P/N 
junction  under  forward  and  reverse 
biased  condition.  Carrier  depletion  in 
reverse  biased  condition. 

Boltzmann  and 
Poisson  equations 

Neutral  ohmic 
region  for  terminal 

Quasi-neutral  doped  region  for  terminal 
with  forward  -biased  P/N  junction 

Diffusion  current 
equation 

High-level- 
injected  ambipolar 
region 

Equal  majority  and  minority  carriers 
which  are  higher  than  the  background 
doping 

Ambipolar  trans- 
port equation 

Current-induced- 

space-charge 

region 

majority  carrier  accumulation  in  the 
region 

Poisson  Equation 

Channel  region 

MOSFET  structure 

MOS  current  equa- 
tion 
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2.3.1  P/N  Junction  Region 


Using  Boltzmann  statistics,  the  carrier  concentration  at 
the  forward  biased  pVn  junction  is  simply  given  as 


Pq  = 


/ 


(2.1) 


where  n,  is  the  silicon  intrinsic  carrier  concentration,  is  the  dop- 
ing concentration  in  the  base,  and  is  the  voltage  drop  across  the 
forward  junction. 

Poisson’s  equation  is  used  to  describe  the  variation  of  junc- 
tion depletion  width  with  voltage 


dV  ^ p(;c) 
dx  e.; 


(2.2) 


where  e,,.  is  the  silicon  dielectric  constant  and  p (r)  is  the  charge  den- 
sity. The  depleted  region  of  the  p"''n  junction  extends  significantly  into 
the  n'  region  during  reverse  bias,  squeezing  the  holes  and  electron 
charges  into  a smaller  volume.  The  positive  charges  in  the  depleted 
region  consist  of  immobile  ionized  donors  and  injected  holes  from  the 
base.  Integrating  (2.2)  twice  over  the  depleted  region  and  rearranging 
the  equation  gives 
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= 


9^epi 


(2.3) 


where  is  the  depletion  region  width,  is  the  built-in  potential  of 
base-collector  junction,  and  ^BC  is  the  voltage  drop  across  the  reverse- 
biased  junction. 


2.3.2  Doping  Region  for  Terminals 

When  the  doping  concentration  of  neutral  ohmic  region  for 
terminal  is  high  (-10^®  cm'^)  as  in  anode  for  the  IGBT,  the  minority 
carrier  (electron)  from  the  epi-base  is  back-injected  to  the  anode.  This 
back-injection  current  is  part  of  anode  current  and  it  can  not  be 
neglected.  The  transport  is  mainly  by  minority  carrier  diffusion  cur- 
rent. The  electron  current  with  negligible  recombination  in  the  p^n 
forward-biased  space-charge  region  is  defined  by  the  recombination  in 
the  quasi-neutral  p"*^  anode  region  as  [ERN91] 

•^BA  ~ JnqPq  (2-4) 


where  = 


^n^A 


electron  diffusion 
in  the  anode. 


, is  the  electron  diffusion  coefficient,  is  the 
length,  and  is  the  acceptor  doping  concentration 


2.3.3  High-Level-Iniected  Ambinolar  Region 


In  general,  the  electron  and  hole  currents  are  given  by 
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(2.5) 


(2.6) 


where  7^,/^  are  the  electron  and  hole  currents,  respectively,  n,p  are 
the  electron  and  hole  carrier  concentrations,  respectively,  A is  the 
device  active  area,  q is  the  electronic  charge,  E is  the  electric  field, 
^„>^p  *ire  the  electron  and  hole  diffusivities,  respectively.  The  first 
terms  in  Eqns  (2.1)  and  (2.2)  are  due  to  drift  and  the  second  terms  are 
due  to  diffusion.  Under  the  high-gain  conditions  of  the  conventional 
bipolar  transistor  analysis,  these  equations  can  be  decoupled  and  the 
transport  of  minority  carriers  in  the  base  for  both  high-  and  low-level 
injection  conditions. 

However,  for  the  low-gain,  high-level  injection  conditions  of 
the  lateral  PNP  BJT/LIGBT  subcell,  the  difference  between  the  elec- 
tron drift  and  diffusion  currents  is  significant  and  the  net  electron 
current  can  not  be  set  to  zero  in  approximating  the  electric  field.  Since 
the  net  electron  current  has  a significant  effect  on  the  hole  drift  cur- 
rent, the  electron  and  hole  transport  equations  can  not  be  decoupled. 
Assuming  quasi-neutrality  and  high-level  injection  in  the  base  region, 
the  current  equations  can  be  written  for  the  ambipolar  transport  case 
by  eliminating  the  electric  field  between  Eqns  (2.5)  and  (2.6)  [HEF86]: 


(2.7) 


(2.8) 
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|X 

where  b = — ia  the  electron/hole  mobility  ratio,  Ij.  is  the  total  current, 
is  the  ambipolar  diffusivity,  which  is  2DnDp/(Dn+Dp).  Equations 
(2.3)  and  (2.4)  depend  on  the  total  current  so  that  the  transport  of  elec- 
trons and  holes  are  coupled. 

The  continuity  equation  for  hole  is 


Xff  q dx 


(2.9) 


where  is  the  high-level  injection  lifetime.  From  Eqn  (2.8)  and  (2.9) 
the  time-dependent  ambipolar  diffusion  equation  is  obtained 


a/  ^ABt 


(2.10) 


where  is  the  ambipolar  diffusion  length  and  is  jD^Xff,  Xjj  is  the 
high-level  injection  lifetime. 

A requirement  in  deriving  this  expression  is  that  the  total 
current  Ij  (anode  current)  is  independent  of  position  in  the  base.  This 
is  satisfied  in  the  LIGBT  because  the  base  current  (electrons)  flows 
from  the  collector  through  the  base  in  the  same  direction  as  the 
injected  hole  current  from  the  emitter.  Equation  (2.10)  does  not  hold 
for  the  traditional  bipolar  transistor  for  which  the  base  current  enters 
from  the  side  in  the  middle  of  the  base. 


2.3.4  Current-Induced-Snace-Charge  Region 

Poisson’s  equation  can  be  used  to  describe  the  voltage  across 


the  CISC  region 
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^ V _ PU)  ^ g_ 


^^EPl 


(2.11) 


where  N^pj  is  the  uniform  doping  concentration  of  the  epi-base  and  the 
n is  the  electron  density  in  the  CISC  region. 

Because  of  relatively  low  doping  in  epi-base  region,  when 
MOS  current  becomes  higher  than  the  critical  current,  7^  = 
the  CISC  region  forms.  Because  the  MOS  current  is  velocity  saturated 
in  the  CISC  region,  the  MOS  current  is  constant  with  x through  the 
entire  CISC  region.  This  is  based  on  the  one-dimensional  high  current 
Kirk  effect  [KIR62]  in  the  bipolar  transistor.  Because  the  Kirk  effect  of 
bipolar  transistor  is  based  on  the  one-dimensional  analysis,  there  can 
be  a error  to  apply  the  one-dimensional  model  into  the  two-dimen- 
sional CISC  region  in  the  LIGBT.  However,  PISCES  simulation  results 
(Fig.  2.5)  show  the  our  one-dimensional  CISC  region  model  can  be 
applied  in  the  LIGBT  CISC  region. 

When  MOS  current  is  less  than  7^,  there  is  no  CISC  region 
near  channel.  Therefore,  equals  to 

From  Fig.  2.5(b),  the  voltage  drop  across  the  CISC  region  as  a 
function  of  x can  be  assumed  to  be  the  grey  triangular  area  based  on 
the  PISCES  simulation,  and  its  expression  is 


= Ie  W . 

CISC  2 


(2.12) 


where  is  the  width  of  the  CISC  region  and  is  assumed  to  be  equal 
to  the  base-collector  depletion  region  width  based  on  the  PISCES  sim- 
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ulation,  is  the  peak  electric  held  in  CISC  region  occurring  at  the 
edge  of  the  LDMOST  channel. 

Poisson’s  equation  for  the  CISC  region,  defining  x=0  as  the 
channel  edge  of  the  base  and  x=Wcisc  as  the  CISC  region  edge  of  the 
quasi-neutral  base,  takes  the  form 


dE  _ q 
~ ^si 


q\ 


MOS 


sat  J 


and 


(2.13) 


f dE  = ^ 
JE  e . 

St 


AT 

<7V  , 

sat  J 


Jo 


From  Eqn.(2.14),  Emax  is 


E = -- 2- 

max  c 
^si 


^epi- 


' MOS 


q\ 


sat  J 


^cisc 


(2.14) 


(2.15) 


Using  Eqns  (2.12)  and  (2.14),  the  final  voltage  drop  across  the 
CISC  region  can  be  expressed  as 


V = -2- 

CISC  Op 
^^si 


'*MOS 
q\ 


sat 


-^epi 


w?... 


where  is  the  current  density  of  the  MOS  current  and  v 
electron  saturation  velocity. 


(2.16) 
is  the 


2.3.5  Channel  Region 

The  MOS  current  of  the  n-channel  LDMOST  is  a nonlinear 
function  of  ^CH  and  the  gate  to  source  voltage  drop  As  the  chan- 
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nels  of  LIGBT  structures  may  be  short  (<2}xm),  current  saturation  due 
to  velocity  saturation  rather  than  vanishing  inversion  charge  and 
channel  length  reduction  during  current  saturation  have  to  be  taken 
into  account.  Because  the  transconductance  parameter  is  different  for 
the  linear  and  the  saturation  region,  both  linear  and  saturation 
transconductance  parameters  must  be  extracted. 

The  gate-to-source  capacitance  and  the  gate-to-drain 

capacitance  c,.  are  both  nonlinear  functions  of  and  On  the 
one  hand,  they  show  the  typical  MOS  capacitance  behavior;  on  the 
other  hand,  charge  sharing  and  charge  shielding  phenomena  giving 
rise  to  a further  degree  of  complexity  and  modeling  effort  have  to  be 
taken  into  account  in  order  to  describe  these  capacitances  accurately. 
The  general  expressions  which  describe  the  current  and  the  charges 
are  presented  in  Section  3.2.2  and  Section  3.3.1. 


2.4  Summary 

The  physical  insights  with  regard  to  two  structures  of  IGBTs 
(LIGBT  and  SOI  LIGBT)  have  been  obtained  by  using  the  two-dimen- 
sional device  simulator  PISCES.  These  insights  were  used  to  identify 
five  different  basic  semiconductor  regions.  The  application  of  correct 
semiconductor  equations  to  each  basic  semiconductor  region  is  critical 
in  modeling  the  power  devices. 

From  the  PISCES  simulations,  we  found  that  IGBTs  have  five 
different  basic  semiconductor  regions:  the  p/n  junction  region,  the  end- 
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doping  region  for  the  terminals,  the  high-level-injected  ambipolar 
region,  the  current-induced-space-charge  region,  and  the  channel 
region.  The  high-level  injection  region  inside  the  IGBT  reduces  the  on- 
state  device  voltage  drop  but  also  increases  the  turn-off  transient.  The 
high-level  injection  region  is  the  area  of  carrier  storage  that  is  a dis- 
tributed phenomenon,  described  by  the  ambipolar  transport  equation. 
Therefore,  the  high-level-injected  ambipolar  region  is  one  of  the  most 
important  regions  in  modeling  the  IGBT.  The  modeling  of  power 
devices  with  a high-level-injection  ambipolar  region  can  be  greatly 
simplified  using  our  modeling  technique.  The  generalized  methodology 
for  the  physics-based  modeling  of  IGBT  developed  herein  can  be  useful 
in  modeling  power  devices  with  wide  and  lightly  doped  regions  such  as 
the  power  bipolar  transistor  and  the  PIN  diode. 


CHAPTER  3 

CELL  MODEL  DEVELOPMENT  FOR  LIGHT 
3.1  Introduction 

This  chapter  presents  a physical  model  for  the  Lateral  Insu- 
lated Gate  Bipolar  Transistor  (LIGHT).  The  LIGHT  is  a high-voltage 
device  that  can  be  integrated  through  both  junction  and  dielectric  iso- 
lation techniques.  The  LIGHT  is  a MOS-controlled  power  switching 
device  for  power  integrated  circuit  applications.  The  schematic  struc- 
ture of  the  device  is  shown  in  Fig.  1.1.  The  n'epi  is  lightly  doped  base 
region  to  support  high  voltages.  The  p'^  (~10^®  cm'^)  cathode  region  has 
a higher  doping  density  than  the  p-body  (~10^®  cm'^)  region  and  is 
included  mainly  to  improve  static  latch-up.  Fig.3.1  also  shows  the 
LIGHT  is  dielectrically  isolated  in  a silicon  tub  from  the  substrate. 
Thus  we  have  not  modeled  the  parasitic  PNP  transistor  formed 
between  the  anode  and  the  p-type  substrate  in  the  conventional 
device. 

The  operation  of  the  LIGHT  simply  can  be  treated  as  a parti- 
tioning of  an  n-channel  LDMOST  (LDMOST/LIGBT  subcell)  and  a 
lateral  PNP  BJT  (lateral  PNP  BJT/LIGBT  subcell).  The  basic  symbolic 
circuit  of  the  LIGHT  is  shown  in  Fig. 3. 1(a).  The  LIGHT  functions  as  a 
bipolar  transistor  that  is  supplied  base  current  by  a MOSFET  [Yil84]. 
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Fig.3.1  (a)  The  symbolic  circuit  of  the  LIGBT. 

(b)  The  basic  equivalent  circuit  of  the  LIGBT. 
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Contrasted  to  VLSI  bipolar  transistors,  the  lateral  PNP  BJT/LIGBT 
subcell  has  a low  base  lifetime,  a low  current  gain  and  is  operated 
mainly  in  high-level  injection  in  the  base  region.  Since  high-level 
injection  prevails  ambipolar  theory  must  be  used  to  describe  the  trans- 
port of  electrons  and  holes  in  the  n'  epi  base  region.  The  modeling  of 
LDMOST/LIGBT  suhcell  is  also  important  because  it  provides  the 
input  base  current  to  the  lateral  PNP  BJT/LIGBT  subcell. 

Although  previous  works  have  modeled  the  operation  of 
Lights  [FOS88,  PAT86,  ERN91],  none  has  specifically  modeled  the  cur- 
rent-induced-space-charge (CISS)  region  in  the  n*  epi  base  region  near 
the  channel.  The  CISC  region  has  a cone-shaped,  cross-sectional  area 
which  is  incorporated  into  the  base  current.  The  partitioning  of  the 
lateral  PNP  BJT/LIGBT  subcell  anode  current  into  a lateral  and  a 
quasi-vertical  part  is  also  important  and  has  not  been  modeled  before. 

We  have  improved  both  subcell  models  by  considering  the 
above  two  criteria.  Modeling  the  CISC  region  in  the  relatively  low- 
doped  base  region  results  in  a more  accurate  LDMOST  model.  Analyti- 
cal geometrical  factors  [CH092]  based  on  the  ambipolar  transport 
equation  for  the  lateral  and  the  quasi-vertical  anode  currents  for  the 
lateral  PNP  BJT/LIGBT  subcell  have  been  derived.  The  total  anode 
current  is  computed  from  two  one-dimensional  currents.  For  the  quasi- 
vertical anode  component,  a geometrical  factor  corrects  for  its  two 
dimensionality. 

The  model  is  implemented  in  the  Saber  circuit  simulator 
[SAB92].  Saber  is  a very  powerful  modeling  tool  for  power  devices 
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which  have  widely  varying  cell  structures  even  within  one  category  of 
device  (for  example,  LDMOST).  The  numerical  “machine*  in  Saber  is 
transparent  to  the  user.  One  template  containing  the  model  equations 
is  all  the  software  that  needs  implementation.  This  is  not  true  of 
SPICE  typically  where,  aside  from  the  model  subroutine,  17  other  sub- 
routines have  to  be  modified. 


3.2  LIGBT  Steady-State  Model 

In  the  on-state  under  steady-state  conduction,  the  LIGBT 
cathode  electrode  is  grounded  and  the  anode  electrode  is  connected 
through  a load  to  a positive  bias  supply.  When  the  gate  voltage  is 
higher  than  the  threshold  voltage,  the  LDMOST  channel  is  inverted. 
The  LDMOST  current  constitutes  the  base  current  of  the  lateral  PNP 
BJT/LIGBT  subcell.  As  the  emitter-base  junction  turns  on,  it  injects 
holes  into  the  epi  base  region  and  thus  modulates  its  conductivity.  The 
injected  hole  current  can  be  divided  into  two  components.  The  first 
component  is  a recombination  current  with  LDMOST  current  injected 
in  the  epi  base  and  anode  regions.  The  second  component  is  collected 
by  the  p'*’  cathode  that  forms  the  collector  of  the  lateral  PNP  BJT/ 
LIGBT  subcell.  The  steady-state  operation  of  the  Lights  can  be  gener- 
ally described  with  a simple  equivalent  circuit,  which  is  shown  in 
Fig. 3. 1(b)  [SIM85,  PAT86,  KIM88].  The  two-dimensional  model  of 
anode  current  for  the  lateral  PNP  BJT/LIGBT  subcell  and  the  model  of 
the  CISC  region  for  the  LDMOST/LIGBT  subcell  are  added  to  the  sim- 
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pie  model  to  improve  the  accuracy  of  the  circuit  simulator  for  the 
LIGHT  devices. 

A system  of  parametric  equations  is  derived  for  the  steady- 
state  electron-  and  hole-current  densities  and  the  excess  carrier  con- 
centration and  the  emitter-base  voltage  in  next  subsection.  These 
equations  are  obtained  by  solving  the  ambipolar  transport  equations 
for  the  boundary  conditions  of  the  lateral  PNP  BJT/LIGBT  subcell. 
The  characteristics  of  the  lateral  PNP  BJT/LIGBT  subcell  are  then 
combined  with  the  LDMOST/LIGBT  subcell  model  to  completely 
describe  the  steady-state  current-voltage  characteristics  of  the  LIGHT. 

3.2.1  Lateral  PNP  BJT/LIGBT  Subcell  Model 

The  analysis  is  performed  using  the  coordinate  system 
defined  in  Fig. 3. 2.  Defining  x = 0 as  the  emitter  edge  of  the  base  and 
X = W as  the  collector  edge  of  the  quasi-neutral  base,  the  steady-state 
boundary  conditions  for  the  hole  distribution  are 

p(W)  =0,  p(0)  = Pq 

Solving  the  steady-state  ambipolar  diffusion  (Eqn(2.10)  with 
^ = 0)  in  the  base  region  with  the  above  boundary  conditions  yields 
one-dimensional  hole  distribution  [FOS86] 


pU)  =Pq 


(3.1) 


36 


X 


0 


W Wb 


Fig. 3. 2 Coordinate  system  used  in  developing  the  LIGHT  model. 
The  junction  space-charge  region  is  indicated  by  the 
dashed  lines. 
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where  Pq  is  the  hole  density  at  the  edge  of  the  emitter-base  junction 
space-charge  region. 

From  the  lateral  PNP  BJT/LIGBT  subcell,  the  anode-cathode 
voltage  (V^jf)  can  be  divided  into  three  components:  (1)  a p^n'  for- 
ward-biased emitter-base  junction  voltage  (V^^),  (2)  a nVp  reverse- 
biased  base-collector  junction  voltage  and  (3)  the  voltage 

drop  across  the  epi  base  region.  Thus, 

^AK  ~ ^EB  ^BC  ^EPl  (3-2) 


Using  quasi-equilibrium  assumption  [MAC87],  Pq  can  be 
expressed  in  terms  of 


(3.3) 


where  is  the  intrinsic  carrier  concentration,  K is  the  Boltzmann 
constant,  and  T is  the  temperature.  W can  be  calculated  through  the 
depletion  approximation  at  the  reverse-biased  collector-base  junction: 


W = W 


B 


qN  . 

^ epi 


(3.4) 


where  Wg  is  the  metallurgical  base  width,  is  the  built-in  potential 
of  base-collector  junction,  is  the  uniform  doping  density  in  the 
base  region,  and  is  the  dielectric  constant  of  silicon. The  base-collec- 
tor reverse-biased  voltage,  will  shrink  the  width  W of  the  quasi- 
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neutral  base.  In  addition,  however,  W is  affected  by  the  current  den- 
sity, which  flows  through  the  depletion  layer  between  base  and 
collector.  This  phenomenon  is  called  Kirk  effect  [KIR62]  in  the  theory 
of  the  bipolar  transistor.  In  the  LIGBT  structure,  the  base-collector 
depletion  layer  spreads  predominantly  into  the  lightly  doped  base,  and 
the  high  collector  current  tends  to  increase  the  width  of  the  quasi-neu- 
tral base  region.  Consequently,  the  width  of  quasi-neutral  base 
becomes 


where  Jj^  is  the  collector  current  density  of  the  lateral  PNP  BJT/ 
LIGBT  subcell,  and  7^  = critical  current,  where  denotes  the 

hole  saturation  velocity.  As  clear  from  Eqn.(3.5),  a high  collector  cur- 
rent will  increase  the  width  of  the  quasi-neutral  base,  and  thus  the 
Kirk  effect  opposes  the  formation  of  the  depletion  region  occurring 
because  of  the  reverse  bias  across  the  base-collector  junction. 


(3.5) 


Now  we  can  determine  with  p^n: 


(3.6) 


(3.7) 


and  total  current  density,  7^ , is 
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Combining  Eqn.(3.6)  and  Eqn.(3.7)  gives 

E = ^ ^KTfb-\yp{x) 

V,,, = i;-£dx 


(3.8) 


(3.9) 

(3.10) 


where  is  defined  at  the  point  where  p(x)  equals 
Finally  V^pj  is 


^EPI  ~ 


y^L^sinh^^^ 


f W-W 

~1T. 


+ 


KTjb-l) 
q ib+l) 


(3.11) 


The  first  of  these  terms  is  due  to  an  ohmic  drop  in  the  epi 
base  region,  whereas  the  second  is  caused  by  the  asymmetric  concen- 
tration gradient  produced  by  the  unequal  electron  and  hole  mobilities 


[GHA77]. 


The  anode  current  is  a sum  of  the  hole  current,  which  is 
injected  from  emitter  (anode)  into  the  epi  base,  and  the  electron  cur- 
rent, which  is  back  injected  into  the  emitter  from  the  epi  base.  For 
negligible  recombination  in  the  emitter-base  junction  space-charge 
region,  can  be  related  to  po  following  conventional  p-i-n-  diode  the- 
ory [GHA77] 


~ JnoPo 


dpjx) 

dx 


x = 0 
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(3.12) 


electron  diffusion  length,  and  is  the  acceptor  doping  concentration 
in  the  anode. 

Under  steady-state  on  condition,  LDMOST  current  flows  into 
the  epi  base  region  through  the  LDMOST  channel.  This  LDMOST  cur- 
rent is  equal  to  the  base  current  of  the  lateral  PNP  BJT/LIGBT 
subcell.  The  base  current  is  the  sum  of  the  recombination  current  in 
the  epi-base  region  and  the  electron  current,  which  is  back  injected 
into  the  emitter  from  the  epi  base.  The  latter  contribution  is  the  same 
as  the  first  term  in  Eqn.(3.12) 


The  collector  current  density  Jj^  is  expressed  as  Previ- 

ously, the  cross  sectional  area  A taken  from  layout  has  been  used  to 
relate  the  terminal  currents  and  Ig  (and  Ij^)  to  the  respective  cur- 


LIGBT devices.  However,  it  has  been  reported  [LIN67]  that  lateral 
devices  with  wide  base  such  as  lateral  PNP  transistor  do  not  have  one- 


(3.13) 


rent  densities.  This  is  true  for  one-dimensional  current  analysis  and 
the  one-dimensional  current  analysis  has  been  used  to  model  the 
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dimensional  current  flow  but  rather  have  a two-dimensional  current 
flow.  To  gain  some  insights  about  the  current  flowline  for  the  LIGBT,  a 
2-D  numerical  device  simulator,  PISCES  has  been  used  (see  Fig.2.6). 

The  PISCES  simulation  implies  the  correctness  in  modeling 
the  anode  current  flow  as  two  components,  a purely  lateral  current 
component  (al  in  the  anode  side-wall  region  and  a quasi-vertical  com- 
ponent along  the  curved  trajectory  7^^,. 

The  one-dimensional  anode  current  based  on  the  ambipolar 
transport  equation  can  be  written  as 


Aao  ~ 


+7«oPo  Pa^j 


(3.14) 


where,  and  p^Xj  is  the  area  of  the  anode  junction  sidewall,  these  geo- 
metrical parameters  are  graphically  shown  in  Fig.  4.1. 

The  purely  lateral  current  which  flows  parallel  to  the  surface 
can  be  modeled  as  a function  of  the  varying  base  width  due  to  lateral 
diffusion  in  the  anode.  The  geometry  of  the  anode  sidewall  is  shown  in 
Fig.3.3.  The  curved  junctions  on  the  anode  side  considered  to  be  cylin- 
drical with  radius  of  curvature  Xj,  the  junction  depth.  The  base  width 
is  a function  of  y 

W^iy)  =W  + Xj(l-cosQ) 
where  6 = 


(3.15) 
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(a) 


(b) 


Fig. 3. 3 (a)  The  geometry  of  anode  and  cathode  sidewall  of  the 
LIGBT.  The  base  width  is  a function  of  x,  is  a base 
width,  and  Xj  is  the  junction  depth,  (b)  close-up  of  geome- 
try of  anode  sidewall. 
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Since  the  lateral  current  component  is  assumed  purely  hori- 
zontal to  the  surface 


(y) 


^■^iNoPoPa^j 


(3.16) 


where 


WAy) 


"A 


dx  = r cosScoth 
Jo 


^W  + ;c^.(l-cos9)^ 


V 


dQ 


y 


(3.17) 


The  lateral  geometrical  factor  is 


(3.18) 


and  it  only  depends  on  the  base  width  W ^ the  junction  depth  Xj,  and 
the  ambipolar  diffusion  length  . 

The  vertical  anode  current  path  can  be  modeled  by  concentric 
curves  with  radius  r as  shown  in  Fig.3.5.  This  model  is  valid  if  the 
maximum  base  width  of  the  device,  shown  in  Fig,  4.1,  is  at  least 
twice  the  length  of  the  device  thickness  Dyj.  as  in  our  case.  However,  if 
Dyj  is  approximately  equal  to  there  is  at  least  a 30%  error  in  the 
vertical  geometrical  factor.  Most  of  conventional  LIGHT  power  device 
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P 


Fig. 3.4  The  geometry  of  the  concentric-curve  model.  The  half-cir- 
cle lines  indicate  the  quasi-vertical  component  of  anode 
current. 
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have  thereby  validating  this  concentric-curve  model  for  the 

LIGBT. 

The  concentric  path  of  quasi-vertical  anode  current  can  be 
seen  and  verified  from  two-dimensional  PISCES  simulation  as  shown 
in  Fig.  2.6.  To  model  the  two-dimensional  anode  current  effects  in  the 
lateral  PNP  BJT/LIGBT  subcell,  the  derivation  of  the  quasi-vertical 
geometrical  factor  is  needed.  Without  this  geometrical  factor,  the 
anode  current  should  be  assumed  to  flow  laterally  through  the  entire 
base  region  in  the  LIGBT.  This  assumption  may  result  the  overestima- 
tion of  the  anode  current  because  more  current  density  can  be  found  in 
the  lateral  current  flow  region  (because  of  short  distance  between  the 
anode  and  the  cathode).  Figure  3.13  shows  the  comparison  of  transient 
simulation  results  between  our  model  and  the  model  without  the  two- 
dimensional  current  effect. 

From  Fig. 3. 4,  the  length  of  the  curvature  is  2r0.  Then  the 
quasi-vertical  geometrical  factor  is 


\\/MAXR  \ 

279'*' 

and 


(3.19) 


COS0  = 


W 


ept 


(3.20) 


where  MAXR  is  the  maximum  radius  and  MINE  is  the  minimum 
radius  of  the  curvature.  Major  parameters,  taken  from  the  layout  and 
cross-sectional  geometry,  are  the  the  epitaxial  thickness 
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and  the  junction  depthoc^..  If  we  assume  the  circle  of  the  MAXR  is  cross- 
ing three  points  as  shown  in  Fig. 3. 4,  then  we  can  find  the  center  point 
of  the  circle.  Finally,  MAXR  is 


MAXR  = 


(3.21) 


and  MINR  is 

MINR  = MAXR-W  . 

ept 


(3.22) 


Then,  Eqn.(3.19)  becomes 


F = ^ f* 

2ac^.Jfl0(l -COS0) 


dQ 


= ^[20  \ + + 


1080  7600 


(3.23) 


where  a = acos 


1- 


J£i 


MAXR-W 


and  b = 


f,  MAXR\ 
acos  1 — — — . 


^epi  )' 


The  total  geometrical  factor  ^GA  is  the  sum  of  and  F^^y. 
Assuming  no  variation  in  V^g  along  the  emitter-base  junction,  the 
complete  analytical  expression  for  the  anode  current  is 

~ ^AO^GA 

where  + 


(3.24) 
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The  ratio  of  and  pQ^y  will  be  used  later  to  calculate  the 
high-level-injection  quasi-vertical  base  region  charges  for  transient 
analysis. 


3.2.2  LDMOST/LIGBT  Subcell  Model 

When  the  gate  voltage  is  raised  above  the  threshold  voltage  of 
the  LDMOST  gate,  the  LDMOST  channel  region  is  inverted.  As  the 
anode  voltage  increases,  the  n-channel  LDMOST  is  turned-on  and  elec- 
tron current  starts  to  flow  into  the  base  region  through  LDMOST 
channel.  The  LDMOST  current  is  the  input  base  current  of  the  lateral 
PNP  BJT/LIGBT  subcell. 

In  the  LIGBT  normal  mode  of  operation,  when  the  LDMOST 
current  is  higher  than  the  critical  current,  7^  = (a  typical 

LIGBT  operating  condition  because  of  the  relatively  low  doping  in  n‘ 
epi  base  region),  the  current-induced-space  charge  (CISC)  region  forms 
in  the  epi-base  region  near  the  channel.  The  voltage  across  the  CISC 
region  is 


V = -F  w 
''CISC  2 


(3.25) 


where  is  the  width  of  the  CISC  region  and  is  assumed  to  be  equal 
the  base-collector  depletion  region  width,  is  the  peak  electric  field 
in  CISC  region  occurring  at  the  edge  of  the  LDMOST  channel.  Then 
becomes 

nxax 
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(3.26) 


The  final  form  of  becomes 


CISC 


CISC 


(3.27) 


As  shown  in  Fig. 2. 5(a),  the  drain  is  located  in  the  base  region 
near  the  channel  and  drain-source  voltage  in  LDMOST/LIGBT 
subcell  is  approximately  in  the  lateral  PNP  BJT/LIGBT  subcell, 
and  the  channel  voltage  is  different  from  The  drain-source 
voltage  is  a sum  of  the  channel  voltage  and  the  voltage  drop  across  the 
CISC  region. 

^CH=  ^Ds-'^cisc  (3.28) 

Equation  (3.26)  explains  why  not  which  is  typically 
used,  must  be  used  to  calculate  an  accurate  value  of  LDMOST  current. 

We  used  the  first  order  MOSFET  current  equation  in  our 
LDMOST/LIGBT  subcell  model  for  its  simplicity.  It  is  suggested  that 
more  refined  LDMOST  current  equations  may  be  needed  in  the 
LDMOST/LIGBT  subcell  model.  When  the  LIGBT  is  in  the  on-state, 
the  LDMOST  is  in  its  linear  region  and  the  MOSFET  first  order  linear 
region  equation  is  given  by 


(3.29) 
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where  transconductance  parameter  for  linear  region,  is  equal  to 

the  product  of  the  oxide  capacitance,  the  surface  mobility,  and  the 
effective  width-to-length  ratio  of  the  LDMOST/LIGBT  subcell.  One  of 
the  unique  properties  of  the  IGBT  is  that  the  transconductance  param- 
eter is  different  for  the  linear  and  the  saturation  region. 

The  LIGBT  is  in  its  current  saturation  region  when  the 
LDMOST  is  in  its  saturation  region.  The  first  order  saturation  current 
equation  of  the  MOSFET  is  given  by 


sat  ^Psat 


'MOS 


( ^GS  ~ ^Th) 


TH> 


(3.30) 


where  Kpsat  is  the  transconductance  parameter  for  saturation  region, 
and  its  extraction  method  will  be  illustrated  in  Chapter  4.  The  satura- 
tion is  assumed  to  be  due  to  velocity  saturation  at  the  end  of  the 
channel  for  simplicity. 

Figure  3.5  shows  the  measured  and  simulated  current-voltage 
characteristics  with  a gate  voltage  of  8 V.  In  the  figure,  the  vertical 
axis  shows  the  anode  current  and  the  horizontal  axis  shows  the  anode- 
cathode  voltage.  Also,  the  I-V  characteristics  for  simulations  of  the 
new  model  using  ^CH  and  the  old  model  using  are  compared. 
Clearly,  the  new  model  provides  a better  fit  to  the  measurement  than 
the  old  model.  The  new  model  shows  a saturation  in  which  is  sup- 
ported by  PISCES  simulations.  Above  = 1.0  V,  the  old  model 
severely  overestimates  the  anode  current  because  continues  to 
increase  monotonically  as  anode-cathode  voltage  increases. 


Ia(A) 
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Fig.3.5  I-V  characteristics  comparison  between  new  model  using 
^CH  model  [FOS88]  using  contrasted  with 

measurement  at  = 8.0. 
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— LIGBT  Steady-State  Model-Routine 

The  flowchart  for  semi-numerical  LIGBT  steady-state  model- 
routine  is  given  in  Fig.  3.6.  The  model  has  three  internal  nodes:  Vm, 
^BC*  ^CH’  solve  for  the  internal  node  voltages,  three  implicit  bi- 
secant semi-numerical  loops  have  been  implemented  in  FORTRAN 
code. 

The  model  routine,  given  flrst  determines  the  node 

voltages  of  the  lateral  PNP  BJT/LIGBT  subcell  (V^^,  Vg^)  using  the 
quasi-static  equations  developed  in  Section  3.2.1.  Vg^  is  assumed  to  be 
^Ds  LDMOST/LIGBT  subcell;  can  be  determined  by  itera- 

tion until  the  LDMOST  current  equals  the  base  current  of  the  lateral 
PNP  BJT/LIGBT  subcell.  Finally,  we  can  determine  Vg.g  using  the  bi- 
secant method. 


3.3  LIGBT  Transient  Model 

In  this  section,  the  turn-off  process  in  the  LIGBT  with  a resis- 
tive load  is  discussed  in  reference  to  the  switching  mode  of  operation. 
After  removal  of  the  gate  voltage,  the  LDMOST  current  falls  rapidly, 
removing  the  base  current  to  the  lateral  PNP  BJT/LIGBT  subcell. 
However,  the  collector  current  of  the  lateral  PNP  BJT/LIGBT  subcell 
falls  more  slowly  because  of  the  predominance  of  the  excess  carrier 
recombination  in  epi  base  region.  Figure  3.7  shows  a circuit  configura- 
tion used  to  measure  the  dynamic  characteristics.  The  measurement  of 
anode  current  turn-off  waveform  of  the  LIGBT  is  shown  in  Fig.  3.8. 
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Fig. 3. 6 The  LIGBT  steady-state  model  flowchart. 
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A6302  Current  Probe 


Fig. 3. 7 Circuit  coniigumtion  used  to  measure  the  dynamic 
characteristics. 


Ia(A) 
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Fig. 3. 8 Measured  LIGBT  anode  current  versus  time  for  turn-off 
transient. 
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The  current  decay  waveform  has  three  phases.  The  sudden  drop  that 
results  from  cessation  of  the  LDMOST  electron  current  is  the  first 
phase.  Following  that,  the  current  decay  is  governed  by  carrier  recom- 
bination in  the  conductivity-modulated  lightly-doped  epi-base  region. 
The  decay  is  slow  because  the  ambipolar  high-level-injection  lifetime 
is  long.  As  the  level  of  injected  carriers  decrease  below  the  doping 
level,  the  current  decay  then  is  dictated  by  the  low-level-injection  life- 
time which  is  shorter  than  that  of  the  high-injection  case. 


3.3.1  Lateral  PNP  B JT/UGBT  Subcell  Model 


Once  the  LDMOST  current  is  removed,  the  stored  excess 
majority  carriers  (electron)  in  the  base  region  decay  by  recombination 
and  by  injection  into  the  anode.  The  excess  of  minority  carriers  (hole) 
in  the  base  continues  to  be  supplied  by  injection  from  anode  and  is 
depleted  by  collection  at  the  base-collector  junction  and  recombination 


in  the  base  in  such  a way  that  quasi-neutrality  is  maintained  through- 
out the  base.  The  collected  hole  current  is  equal  to  the  lateral  PNP 
BJT/LIGBT  subcell  current  during  the  slowly  decaying  current  phase 
at  a constant  voltage  shown  in  Fig.  3.8. 

The  predominant  charges  stored  in  the  lateral  PNP  BJT/ 
LIGBT  subcell  are  Qbe>Qbc>  »nd  Qj^.  Q^e  and  Q^c  are  the  hole 
charges  stored  in  the  quasi-neutral  base  of  the  lateral  PNP  BJT/ 
LIGBT  subcell  associated  with  the  anode  and  cathode,  respectively.  In 
general,  and  Q^c  depend  on  both  the  and  This  means 


dQ 


that  the  charging  current  has  two  terms 


dt 


56 


dV^c) 


(3.31) 


The  first  term  can  be  represented  as  a capacitor,  and  the  sec- 


transient  model,  the  simulation  of  turn-off  transient  will  be  in  signifi- 


time  in  the  quasi-neutral  base. 

Since  the  injected  base  region  charge  is  a function  of  anode 
current,  each  and  should  be  partitioned  into  two  components, 
lateral,  and  quasi-vertical.  The  partitioning  is  defined  by  integrating 
the  carrier  continuity  equations  in  the  base.  This  is  physically  repre- 
sentative of  the  bipolar  transistor  [FOS86] 


where  and  lateral  components  of  the  base  charge 

associated  with  anode  and  cathode,  respectively,  Xj  is  the  junction 
depth,  and  is  the  device  width,  shown  in  Fig.  4.1. 


ond  term  as  a transcapacitance.  Without  the 


term  in  the  LIGHT 


(3.32) 


(3.33) 


From  Eqn.  (3.1) 


p(x,t)  =p(0,0 


(3.34) 
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Inserting  Eqn.  (3.32)  into  Eqn.  (3.30)  and  Eqn.  (3.31)  yields 
Qbel  = ^^vw^j^aP  (0.  t)  j (3.35) 

and 

Qbcl  ~ (3.36) 


The  time  dependence  of  and  are  defined  by  the 

time  dependences  of  the  emitter-base  and  base-collector  voltages 
which  control  p(0,t)  and  W. 

The  geometrical  factors  can  be  used  to  calculate  Qg^y  and 
Qbcv^  quasi-vertical  components  of  the  base  charge  associated  with 
and  cathode,  respectively 


O = 

^BEV  - p '^BEL 
^ GAL 


O = 

- p ^BCL 
^ GAL 


(3.37) 

(3.38) 


Since  and  are  the  sum  of  its  lateral  and  quasi-verti- 
cal components 

^BE  ~ Qbel'^Qbev  (3.39) 

Qbc  ~ ^BCL’^O-Bcy  (3.40) 

Qjq  is  the  reverse-biased  base-collector  depletion  region 


charge 
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Qjc  = 


(3.41) 


where  A is  the  cross-sectional  area,  is  the  base-collector  depletion 
region  width. 


LDMQST/LIGBT  Suhcell  MoHpI 

The  charges  in  the  LDMOST/LIGBT  subcell  are 

*'Cro  ^ ^GD  ^ 

and  Qc/sc-  LDMOST  charge  model  we  did  not  include  the  chan- 

nel charge  in  and  for  simplicity.  It  is  suggested  that  our 
LDMOST  charge  model  should  be  refined  to  improve  the  accuracy  of 
the  transient  analysis. 

Qgs^  the  gate  oxide  capacitance  of  the  source  overlap  charge 

^GS  ~ ^soPvW^ox^GS  (3.42) 

where  is  the  length  of  the  gate  oxide  overlapping  in  source. 

Qgd  gate-drain  capacitance  times  gate-drain  voltage 

^GD  ~ ^GD^GD  (3.43) 

Cgd  is  equal  to  the  gate  oxide  capacitance  of  the  gate-drain 
overlap  for  y^s^^GS*  ^ds^^gs  silicon  beneath  the  gate- 

drain  overlap  becomes  depleted  and  the  gate-drain  capacitance  con- 
sists of  the  series  combination  of  the  gate-drain  overlap  oxide 
capacitance  and  the  gate-drain  overlap  depletion  capacitance. 


C 


GD 


^oxd 
^oxd^gdj 
^^oxd  ^gdp 


for  ( ^ Vq^) 
for  { Vpg  > Vq^) 


(3.44) 
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where  - ^^oiDywCox  is  gate  oxide  capacitance  of  the  drain  overlap 
charge  and  i®  the  gate-drain  depletion  capaci- 

tance, is  the  length  of  the  gate  oxide  overlapping  in  drain, 
the  gate-drain  depletion  width,  is  proportional  to  the  square-root  of 
the  gate-drain  voltage. 

Unlike  the  VIGBT  [HEF90],  the  silicon  beneath  the  gate- 
drain  overlap  not  only  becomes  depleted  but  also  becomes  a excess 
electron  accumulation  region  in  the  LIGBT.  The  accumulation  region 
is  called  the  current-induced-space-charge  (CISC)  region.  After  the 
gate  is  turned-off,  the  charges  in  CISC  region  dissipate  by  recombina- 
tion in  the  base  region.  The  charge  in  the  CISC  region  is 


^CISC  ~ 9^cisc^VW\ 


i:SA5) 


is  the  cross-sectional  area  of  CISC  region. 


a.3.3  LIGBT  Network  Representation 

The  complete  network  representation  of  the  new  LIGBT 
model  is  shown  in  Fig.  3.9.  The  model  is  not  a simple  equivalent  cir- 
cuit. The  elements  shown  represent  physical  components  of  current 
and  voltage  in  the  LDMOST,  described  implicitly  by  the  analysis, 
which  requires  a numerical  solution. 

The  transient  behavior  of  the  LIGBT  is  simulated  using  both 
the  quasi-static  currents/voltages  and  charging  currents  The 
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G 


Fig.3.9  Complete  network  representation  of  the  charge-based 
LIGBT  model. 
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charges  within  the  device  are  expressed  in  terms  of  node  voltages,  and 
the  transient  charging  currents  are  represented  by  time  derivatives  of 
the  charges. 


3.4  LIGBT  Model  Verification 

Measured  and  simulated  current-voltage  characteristics  for 
three  different  values  of  gate  voltages  are  shown  in  Fig.  3.10.  The 
measured  data  in  Fig. 3. 10  were  taken  with  an  HP-4145  parameter 
analyzer.  The  following  features  of  the  static  characteristics  are  indi- 
cated in  Fig.  3.10:  1)  the  diode  voltage  offset  due  to  the  anode-epitaxial 
layer  of  p/n  junction  at  low  V^k  dess  than  -0.7V),  2)  the  on-state 
region  which  has  low-resistance  due  to  conductivity  modulation,  3)  the 
current  saturation  region  due  to  saturation  of  the  LDMOST  channel 
current,  and  4)  the  triode-like  region  which  is  referred  to  as  the  linear 
region  for  LDMOST.  The  agreement  between  measurement  and  simu- 
lation for  1)  the  diode  voltage  offset  validates  the  implementation  of 
Veb,  2)  the  on-state  characteristics  validates  the  implementation  of 
Vepi,  3)  the  saturation  characteristics  validates  the  implementation  of 
^MOS»  Ik>  1^,  and  4)  the  triode-like  region  validates  the  implemen- 
tation of  Vcfj.  The  agreement  is  good  even  without  the  optimization  of 
the  parameters. 

Figure  3.11  compares  the  channel  voltages  from  PISCES 
simulation  and  our  model  as  a function  of  , contrasted  with  the 
characteristic  for  ^DS  • Because  the  value  of  is  assumed  to  be 


Ia,  amps 
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Fig. 3. 10  Measured  and  simulated  I-V  characteristics  for  four  dif- 
ferent values  of  gate-source  voltages  (Vqs=6.0,  8.0,  10.0, 
and  12.0  volts) 


VOLTAGE  (V) 
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Fig. 3. 11  PISCES  and  model  simulated  versus 

trusted  with  at  = lO.Ovolts . 


con- 
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equal  to  the  value  of  increases  monotonically  as 

increases.  We  used  in  our  model  and  obtained  poor  agreement 
with  the  measurement  data  in  Fig.  3.9,  which  leads  us  to  believe  the 
model  optimization  is  suspected  in  earlier  models  [FOS88,  PAT86, 
ERN91].  We  emphasize  we  did  not  perform  any  optimization  of  the 
model  parameters  for  the  simulation  in  Fig.  3.9.  Both  PISCES  and  our 
model  simulations  show  the  saturation  of  as  increases.  It  is 
obvious  that  a more  accurate  value  of  the  LDMOST  current  can  be 
obtained  by  using  instead  of  . Since  the  LDMOST  current  is 
the  input  base  current  of  the  lateral  PNP  BJT/LIGBT  subcell  and  of 
comparable  magnitude  to  anode  current,  an  accurate  value  for  the 
LDMOST  current  is  essential  in  predicting  an  accurate  value  for  anode 
current. 

The  two-dimensional  model  of  anode  current  also  improves 
the  accuracy  of  our  LIGBT  model.  The  lateral  and  quasi-vertical  geo- 
metrical factors  vary  dependent  on  bias  conditions.  At  = 3 V and 

^Gs  ~ ^GAL  ^GAV  0.84.  In  the  conventional  one- 

dimensional model,  the  anode  current  was  calculated  simply  as  the 
anode  current  density  times  area.  Thus,  it  would  require  optimization 
in  the  parameters  for  the  conventional  one-dimensional  model  to  lit 
the  simulation  result  to  the  measurement. 

In  Fig.  3.12,  a constant-anode-supply-voltage  test  circuit  is 
used  to  verify  the  transient  model.  Measured  and  simulated  transient 
turn-off  characteristics  of  the  LIGBT  are  compared  in  Fig.  3.13.  The 
turn-off  transient  simulation  shows  good  agreement  with  measure- 
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Fig. 3. 12  Circuit  configuration  of  an  LIGBT  with  a constant  anode 
voltage. 


Anode  Current  (Amps) 
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Fig. 3. 13  Measured  and  simulated  LIGHT  anode  current  versus 
time  for  turn-off  transient  at  Vgs=10  V and  Vak=  2.0  V. 
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ment.  The  high-level  injection  lifetime  of  5.6^sec,  the  threshold 
voltage  of  3.4volts,  and  device  width  of  120|im  were  used  to  simulate 
the  transient  analysis.  The  high-level  injection  lifetime  was  calculated 
from  the  transient  measurement  (Fig.  3.8)  and  the  parameter  extract- 
ing method  is  well  explained  in  Chapter  4.  The  good  fit  between  the 
simulation  and  the  measurement  on  1)  initial  rapid  drop  validates  the 
implementation  of  LDMOST  current  Imos>  fli^d  2)  the  slow  decay  of  the 
turn-off  current  waveform  validates  the  model  for  and  the  lat- 
eral and  the  quasi-vertical  two-dimensional  base  region  charges,  and 
the  partitioned  charges,  and 

In  Fig.  3.13  we  have  the  simulated  result  (solid  curve)  which 
have  been  obtained  using  our  model.  So  here  the  CISC  region  model 
and  the  two-dimensional  current  model  have  been  included  in  the  cal- 
culations. The  dashed  curve  represent  the  simulated  result,  which  is 
obtained  without  the  CISC  region  model  and  the  two-dimensional  cur- 
rent model,  but  the  other  details  of  the  theory  have  been  unchanged. 

There  are  two  distinctive  differences  between  two  curves. 
First  is  the  total  current  magnitude  difference.  As  stated  in  Section 
3.2.2,  the  anode  current  can  be  overestimated  without  using  the  CISC 
region  model.  Without  the  two-dimensional  current  model  the  anode 
current  also  can  be  overestimated  because  the  anode  current  is 
assumed  to  flow  linearly  in  whole  area  of  the  device.  Second  is  the  dif- 
ference of  turn-off  decay  time.  The  turn-off  decay  time  of  the  solid 
curve  (~7}i,s)  is  much  faster  than  the  turn-off  decay  time  of  the  dashed 
curve  (~12|is).  The  turn-off  time  of  4ps  is  more  reasonable  value  than 


Anode  Current  (Amps) 
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Time  (sec) 


Fig. 3. 14  Turn-off  transient  anode  current  simulation  comparison 
between  our  LIGBT  model  (solid  curve)  and  the  LIGHT 
model  without  the  CISC  region  model  and  the  two-dimen- 
sional current  model  (dashed  curve). 
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the  turn-off  time  of  7)is  for  a high-level  injection  lifetime  of  5.6|xs 
because  of  the  significant  carrier  removal  via  recombination  in  the  p+ 
anode.  The  dashed  curve  has  longer  turn-off  decay  time  due  to  its  one- 
dimensional current  model.  Without  two-dimensional  current  model, 
the  base  charge  is  calculated  from  the  integration  of  hole  charge  from 
0 to  w (base  width)  times  device  area.  This  can  lead  the  overestimation 
of  the  base  charge. 


3.5  Summary 

We  have  developed  an  accurate  steady-state  and  transient 
model  for  LIGBT.  Improvements  on  both  subcell  models  of  the  LIGBT 
have  resulted  in  accurate  simulation  results  without  parameter  opti- 
mization. In  the  LDMOST/LIGBT  subcell,  we  have  modeled  the  CISC 
region  based  on  PISCES  simulation.  The  model  of  CISC  region  has 
been  shown  to  be  of  great  importance  for  improving  the  accuracy  of  the 
LDMOST  model.  In  the  lateral  PNP  BJT/LIGBT  subcell,  we  have 
derived  analytical  geometrical  factors  based  on  the  ambipolar  trans- 
port equation.  The  lateral  and  quasi-vertical  geometrical  factors 
enhance  the  analysis  of  the  two-dimensional  anode  current,  and,  thus, 
result  in  the  accurate  modeling  for  steady-state  and  transient  condi- 
tions for  the  LIGBT.  The  results  of  the  simulation  show  good  fit  with 
both  steady-state  and  transient  measurements  without  the  optimiza- 
tion of  the  model  parameters.  The  LIGBT  model  has  been  used  to 
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simulate  various  external  circuits  and  has  been  proven  valid  for  simu- 
lating various  power  circuits.  This  will  be  shown  in  Chapter  5. 


CHAPTER  4 

PARAMETER  EXTRACTION  FOR  LIGHT 
4.1  Introduction 

Contrasted  to  microelectronic  devices  where  the  steady-state 
current-voltage  characteristics  in  conjunction  with  interelectrode 
capacitance-voltage  characteristics  are  sufficient  to  extract  most  of  the 
device  model  parameters,  the  dynamic  characteristics  must  in  general 
be  examined  to  characterize  the  devices  and  to  extract  model  parame- 
ters for  merged  power  devices  such  as  the  IGBT.  This  is  because  the 
characteristics  of  the  internal  MOSFETs  and  bipolar  transistors  are 
convoluted  in  the  steady-state  characteristics.  The  dynamic  measure- 
ments require  high  precision  and  the  ability  to  make  calculations  on 
waveforms  obtained  from  measurements.  In  addition,  calculations  are 
required  to  separate  the  measured  characteristics  into  the  characteris- 
tics of  the  internal  MOSFETs  and  bipolar  transistor. 

In  our  LIGHT  model,  parameters  can  be  divided  into  two 
major  categories.  The  first  category  is  device  parameters.  The  device 
parameters  are  layout  dependent  and  thus  can  differ  from  device  to 
device  in  the  same  IC.  The  LIGHT  model  has  nine  device  parameters 
defined  in  Table  4.1  and  Fig.  4.1.  The  second  category  is  model  parame 
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Table  4.1  LIGBT  DEVICE  PARAMETERS 


Name 

Description 

Units 

Default 

Dyw 

Device  Width 

m 

240e-6 

Dvt 

Device  Thickness 

m 

18.0e-6 

Xj 

Junction  Depth 

m 

6.5e-6 

Pa 

Anode  Perimeter  Length 

m 

18.0e-6 

Laj 

Anode  Implantation  Win- 
dow Length 

m 

6.0e-6 

WfiM 

Maximiun  Base  Width 

m 

65.0e-6 

Wbo 

Minimum  Base  Width 

m 

50.0e-6 

Lsol 

Gate  Oxide  Source  Over- 
lapping 

m 

l.Oe-6 

^DOL 

Gate  Oxide  Drain  Overlap- 
ping 

m 

4.0e-6 

Table  4.2  LIGBT  MODEL  PARAMETERS 


Name 

Description 

Units 

Default 

High-Injection  Lifetime 

sec 

2.4e-6 

^TH 

Threshold  Voltage 

V 

4.8 

Nepi 

Uniform  Base  Doping 

m-3 

1.0e21 

Adp 

Peak  Anode  Doping 

m-3 

1.0e25 

Pbd 

P-Body  Doping 

m-3 

5.0e22 

Xpsat 

Saturation  Transconduc- 
tance Parameter 

A/W^ 

4.1e-4 

Xpiin 

Linear  Transconductance 
Parameter 

A/W^ 

9.6e-4 
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ters.  The  model  parameters  are  process-dependent  and  generally  com- 
mon to  all  the  LIGBT  devices  on  the  same  chip.  The  complete  set  of 
seven  parameters  is  shown  in  Table  4.2.  The  key  model  parameters  are 
VxH,  Kpsat,  and  Kp^n,  and  can  be  obtained  by  DC  and  turn-off  tran- 
sient measurement.  The  others,  Ngpj,  A^p  and  Pqd,  are  technological, 
and  can  be  obtained  from  processing  information. 

The  PISCES  simulation  results  are  used  to  extract  the 
parameters  because  we  could  not  get  the  good  measurement  data  of 
transient  analysis.  Three  basic  types  of  simulations  are  used  to  extract 
the  key  model  parameters:  1)  The  turn-off  current  tail  decay  rate  ver- 
sus anode  current  characteristic  is  used  to  extract  the  high-injection 
lifetime  Th.  2)  The  saturation  current  versus  gate  voltage  characteris- 
tic is  used  to  extract  the  MOSFET-channel  saturation-region 
transconductance  parameter  Kps^t  and  MOSFET-channel  threshold- 
voltage  VxH-  3)  The  on-state  voltage  versus  gate  voltage  characteristic 
at  a constant  anode  current  is  used  to  extract  the  MOSFET-channel 
linear-region  transconductance  parameter  Kpun-  The  three  simula- 
tions are  performed  in  the  indicated  sequence  because  the  extracted 
parameters  taken  from  simulated  data  require  the  values  of  the  previ- 
ously extracted  parameters. 


4.2  High-Level-Iniection  Lifetime 

The  simulated  turn-off  transient  anode  current  of  the  LIGBT 
with  the  constant  anode  supply  voltage  test  circuit  is  shown  in  Fig. 
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Fig. 4.1  LIGBT  structure  with  device  parameters 
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4.2.  Two  different  phases  of  the  turn-off  current  are  shown.  An  initial 
fast  drop  that  results  from  the  termination  of  LDMOST  channel  cur- 
rent followed  by  slow  decay  due  to  the  recombination  of  excess  holes 
injected  by  the  anode  into  the  base  region.  The  slow  decay  also  has  two 
phases,  high-level-injection  lifetime  decay  and  low-level-injection  life- 
time decay.  During  the  first  phases  of  slow  decay,  the  current  decay  is 
dictated  by  the  high-level  injection  lifetime.  As  the  level  of  injected 
carriers  decreases  below  the  base  doping-level,  the  lifetime  decreases 
because  the  base  enters  the  low-level  injection  condition.  Therefore, 
during  the  second  phase  of  the  slow  decay,  the  current  decay  is  dic- 
tated by  the  low-injection  lifetime  which  is  shorter  than  that  for  the 
high-injection  case  [DAR87]. 

The  first  phase  of  the  slow  decay  has  been  matched  to  an 
exponential  waveform  with  99.998%  of  accuracy.  If  we  assume  the  first 
phase  of  the  slow  decay  has  an  exponential  form,  then  the  anode  cur- 
rent is 


where  Iq  is  the  initial  value  of  anode  current  of  the  first  phase  of  the 
slow  decay. 


(4.1) 


(4.2) 
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Fig, 4. 2 Simulated  turn-off  transient  anode  current  versus  time 
with  the  constant  anode  supply  voltage  test  circuit. 


-[Ln(I(t)/Ic)] 
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Fig. 4. 3 The  anode  current  tail  decay  rate  indicating  that  the  slope 
of  the  least-square  fit  is  equal  to  the  1/xn- 
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Figure  4.2  shows  the  simulated  values  of  the  turn-off  tran- 
sient of  the  anode  current  versus  time,  and  a least-squares  fit  to  the 
anode  current  decay.  The  value  of  high-level  injection  lifetime,  thi  can 
be  extracted  from  the  slope  of  the  least-square  fit  of  Eqn.(4.2).  The 
graphic  output  of  this  extraction  routine  is  given  in  Fig. 4. 3 and  the 
slope  is  equal  to  the  I/th- 


4.3  Threshold  Voltage  and  Saturation  Transconductance 

Figure  4.4  shows  a simulated  constant-voltage  turn-off  tran- 
sient current  waveform  with  four  different  gate-source  voltages  at  a 
constant  anode  voltage.  The  characteristic  consists  of  an  initial  rapid 
fall  followed  by  a slowly  decaying  current  tail.  Because  the  source  of 
the  majority  carriers  (the  LDMOST  current)  is  absent  during  the 
slowly  decaying  phase  of  the  current  waveform,  the  total  stored  excess 
majority  carriers  in  the  base  decays  by  recombination  and  by  injection 
into  the  anode.  The  excess  majority  carriers  in  the  base  continue  to  be 
supplied  by  injection  from  the  anode  and  are  depleted  by  collection  at 
the  base-collector  junction  and  recombination  in  the  base  in  such  a 
way  that  quasi-neutrality  is  maintained  throughout  the  base.  The  col- 
lected minority  carrier  (hole)  current  is  equal  to  the  total  LIGBT 
current  during  the  slowly  decaying  current  phase  at  a constant  volt- 
age. It  is  shown  that  the  initial  rapid  fall  in  current  is  not  equal  to  the 
base  current  (LDMOST  current)  for  two  reasons:  1)  the  anode  voltage 
and  the  quasi-neutral  base  width  are  in  general  different  during  the 
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Fig. 4. 4 Simulated  constant  voltage  turn-off  current  waveform 
with  four  different  gate-source  voltages  at  a constant 
anode  voltage  Vak=2.0. 
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collector  current  decay  than  during  steady-state  [HEF88],  and  2) 
removal  of  the  electron  current  also  reduces  the  hole  current.  This 
indicates  that  the  quasi-static  approximation  is  not  valid  for  the 
devices  governed  by  ambipolar  transport  equation.  Hence,  the  initial 
drop  of  the  forward  current  AI  is  usually  less  than  the  LDMOST  chan- 
nel current  due  to  the  compensation  made  by  the  base-collector 
junction  displacement  current.  However,  the  effect  of  those  compensa- 
tion currents  can  be  minimized  provided  that  we  limit  Vce  to  a 
relatively  small  value  by  selecting  a small  resistor  in  our  measure- 
ment [SHE93].  Therefore,  AI  can  be  approximately  considered  as  Imos- 
The  common-collector  current  gain  P is  given  by 


If-Ic 


= constsal 


(4.3) 


When  the  LDMOST  current  Imos  ©inters  the  current  satura- 
tion regime,  the  base  current  of  the  lateral  PNP  BJT/LIGBT  subcell 
becomes  constant.  Therefore,  the  LIGBT  anode  current  saturates  when 
the  internal  LDMOST  current  saturates  where  the  LDMOST  satura- 
tion current  is  magnified  (1+P)  by  the  steady-state  common  collector 
current  gain  of  the  lateral  PNP  BJT/LIGBT  subcell.  Fig.4.5  shows  the 
lateral  PNP  BJT/LIGBT  subcell  common-collector  current  gain  P ver- 
sus anode  current,  which  was  extracted  from  the  dynamic 
characteristics  of  the  LIGBT.  At  high  currents,  the  value  of  p 
decreases  due  to  the  increased  rate  of  electron-current  injection  into 
the  anode.  Once  P is  known,  the  saturation  characteristics  of  the  inter- 
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Anode  current  (A) 


Fig. 4. 5 The  common  collector  current  gain  versus  anode  current 
for  the  lateral  PNP  BJT/LIGBT  subcell  of  LIGBT. 
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nal  LDMOST  are  obtained  from  LIGBT  saturation  characteristics  by 
dividing  out  the  lateral  PNP  BJT/LIGBT  current-gain 


sat 


f. 


at 

MOS 


(i  + P) 


(4.4) 


The  square  root  of  the  saturation  current  is  linearly  related 
to  the  gate-source  voltage  with  a zero  current  intercept  of  Vth 
with  a slope  of  JKp^^/2  such  that 


(4.5) 


Figure  4.6  shows  the  square  root  of  the  LDMOST  saturation 
current  versus  gate-source  voltage  obtained  by  dividing  the  simulated 
LIGBT  saturation  current  by  the  common-collector  current  gain.  The 
values  of  Kp^^t  and  Vth  are  extracted  from  the  slope  and  the  zero-cur- 
rent intercept  of  this  curve. 


4.4  Linear  Transconductance 

Figure  4.7  shows  the  steady-state  linear  characteristics  of  the 
LIGBT,  indicating  the  on-state  voltage  for  different  gate  voltages  at  a 
constant  anode  current.  From  the  lateral  PNP  BJT/LIGBT  subcell,  the 
anode-cathode  voltage  ( can  be  divided  into  three  components:  1)  a 
pVn'  forward-biased  emitter-base  junction  voltage  (V^g),  2)  a nVp 
reverse-biased  base-collector  junction  voltage  (V^^),  and  3)  the 
voltage  drop  across  the  epi-base  region 


(LDMOST  Saturation  Current) 
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Fig. 4. 6 The  square  root  of  the  LDMOST  saturation  current  indi- 
cating that  the  slope  is  used  to  extract  Kpg^t  and  the  zero 
current  intercept  is  used  to  extract  Vth- 


Anode  Current  (A) 
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Anode  Voltage  (V) 


Fig.4.7  The  steady-state  linear  characteristics  of  the  LIGBT,  indi- 
cating that  the  on-state  voltage  for  different  gate 
voltages  at  a constant  anode  current. 
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^AK  ~ ^EB  ^BC  ^EPI  (^-6) 

For  the  LIGBT  on-state,  the  LDMOST  channel  current  is  in 
the  linear  region,  and  this  LDMOST  channel  current  is  magnified  by 
the  common  collector-current  gain  of  the  lateral  PNP  BJT/LIGBT  sub- 
cell shown  in  Fig.  4.5 

Jin 

Jin  •j 

^MOS  = (I+P)  (4-7) 

The  drain-source  voltage  in  LDMOST/LIGBT  subcell  is 
approximately  in  lateral  PNP  BJT/LIGBT  subcell,  and  if  it  is 
assumed  that  the  LDMOST  channel  current  is  given  by  the  first  term 
in  the  linear  region  expression  of  Eqn.(3.27) 

K.  = V'm  + + (l+P)/C„Jvoj-V,.J  (“'S) 

From  Eqn.(4.8),  it  is  expected  that  the  on-state  voltage  versus 
1/(Vgs"Vth)  a constant  anode  current  is  linear  with  a slope  of  It/ 
^Plin(l+P)  with  an  infinite  gate  voltage  intercept  of  Veb+Vepi- 

Fig.  4.8  shows  the  linear  graph  of  anode  voltage  versus  I/CVqs-Vth)- 
As  indicated  in  Fig. 4. 8,  the  value  of  Kpij^  is  obtained  from  the  slope  of 
the  least-squares  fit  to  the  measured  on-state  voltage  versus  1/(Vgs" 
Vth)  at  a constant  anode  current  over  the  range  of  gate  voltages  where 
the  LDMOST  channel  current  is  given  by  the  first  term  in  the  linear 
region  expression  of  Eqn.(3.27).  The  value  of  Kpji^  is  thus  obtained 
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Fig. 4. 8 On-state  voltage  versus  1/(Vgs-Vth)  at  a constant  anode 
current  over  the  range  of  gate  voltages.  The  slope  is  used 
to  extract  Kpij^. 
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from  tliis  oxtr&ction  stop  using  tlio  v&luos  of  Vfjj  snd  P c&lcul&tod  from 
model  parameters  obtained  in  previous  extraction  steps. 


4.5  Summary 

The  parameter  extraction  method  for  the  LIGBT  from  simula- 
tion results  have  been  developed  and  described  in  this  chapter.  It  has 
been  shown  that  unlike  parameter  extraction  for  microelectronic 
devices,  the  dynamic  characteristics  must  be  used  because  the  charac- 
teristics of  the  internal  LDMOST  and  lateral  PNP  bipolar  transistor 
subcells  are  convoluted  in  the  steady-state  terminal  characteristics  of 
merged  power  devices  such  as  the  LIGBT.  The  dynamic  waveforms  con- 
tain many  features  that  isolate  different  model  parameters.  The  model 
parameters  are  obtained  sequentially  by  selecting  features  of  the 
device  characteristics  that  isolate  parameters,  and  using  the  parame- 
ters obtained  from  previous  extraction  steps  to  calculate  the  model 
parameters  from  the  simulated  characteristics. 


CHAPTER  5 

MODEL  IMPLEMENTATION  AND  VERIFICATION 

5.1  Introduction 

A physics-based  model  for  the  low-gain,  high-level  injection 
lateral  PNP  BJT/LIGBT  subcell  combined  with  the  LDMOST/LIGBT 
subcell  model  has  been  developed,  resulting  in  a general  purpose 
model  for  the  LIGBT.  The  model  has  been  verified  for  various  external 
circuit  conditions  and  for  the  full  range  of  static  and  dynamic  condi- 
tions in  which  the  LIGBT  is  intended  to  be  operated. 

The  purpose  of  this  chapter  is  to  describe  the  methodology  for 
implementing  the  physical  model  of  the  LIGBT  into  the  Saber  circuit 
simulator  [ANA92],  and  to  provide  a Saber  LIGBT  model  that  can  be 
used  for  general  purpose  circuit  simulations.  Furthermore,  the  tech- 
niques demonstrated  for  modeling  the  LIGBT  with  the  Saber  circuit 
simulator  are  generally  applicable  to  modeling  other  power  semicon- 
ductor devices.  Many  diverse  types  of  power  devices  are  currently 
available  with  structures  that  differ  significantly  from  one  another, 
and  different  device  model  equations  are  generally  required  to 
describe  each  device  type.  This  is  in  contrast  to  microelectronic  devices 
for  which  many  different  devices  can  be  described  by  using  the  appro- 
priate model  parameters  in  standard  device  model.  This  is  the  reason 
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that  Saber  became  a very  powerful  modeling  tool  for  power  devices 
which  have  widely  varying  cell  structures  even  within  one  category  of 
device.  The  numerical  machine  in  Saber  is  transparent  to  the  user, 
and  one  template  containing  the  model  equations  is  the  only  imple- 
mentation needed.  This  is  not  true  of  SPICE  typically  where,  aside 
from  the  model  subroutine,  17  other  subroutines  have  to  be  modified. 
The  techniques  described  in  this  chapter  for  implementing  the  LIGBT 
model  equations  into  the  Saber  circuit  simulator  are  essential  for 
accurate  simulation  of  power  semiconductor  devices.  The  model  imple- 
mentation into  the  Saber  circuit  simulator  closely  follows  the 
treatment  given  by  Hefner  and  Diebolt  [HEF91],  we  will  recapitulate 
the  whole  development  for  easy  reference. 


5.2  Implementing  the  LIGBT  Model  into  Saber 

To  describe  the  behavior  of  a system  such  as  an  electrical  net- 
work using  the  Saber  circuit  simulator,  the  interconnections  of  the 
different  components  of  the  system  are  described  using  a network  list- 
ing (net-list).  The  net-list  contains  a statement  for  each  component  of 
the  system  that  defines  the  name  of  the  model  template  used  to 
describe  the  component,  the  terminal  connection  points  of  the  compo- 
nent, and  the  values  of  the  model  parameters  that  are  to  be  changed 
from  the  default  values  of  the  generic  model  template.  The  model  that 
describes  each  of  the  components  of  the  system  can  be  accessed  from 
the  Saber  libraries  of  standard  component  models  or  from  user-defined 
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Saber  templates  where  the  equations  that  describe  the  physical  behav- 
ior of  the  device  are  implemented.  The  implementation  of  the  LIGBT 
model  equations  into  a Saber  template  is  described  in  this  section,  and 
the  net-lists  describing  the  operation  of  the  LIGBT  within  various  test 
circuits  are  given  in  the  next  section. 

5.2.1  Saber  Tempi  ate 

Saber  templates  are  written  in  the  MAST  modeling  language 
which  is  similar  to  the  C programming  language  with  the  addition  of 
specially  designed  modeling  constructs  which  facilitate  the  implemen- 
tation of  Kirchhoff’s  law  and  aid  convergence.  Both  user-defined 
models  and  the  standard  Saber  library  models  are  implemented  in 
Saber  templates  using  the  MAST  modeling  language.  Electrical  compo- 
nent models  are  implemented  into  templates  by  expressing  the  current 
through  each  element  of  the  component  in  terms  of  the  system  vari- 
ables of  the  component:  system  variables  for  electrical  component 
models  consist  of  terminal  node  voltages,  internal  node  voltages,  and 
explicitly  defined  system  variables.  The  simulator  solves  for  the  sys- 
tem variables  of  the  entire  network  such  that  the  net  current  into  each 
node  of  the  system  sums  to  zero,  (i.e.,  Kirchhoff’s  current  law  is  satis- 
fied), and  such  that  the  equations  defining  the  explicitly  defined 
system  variables  for  each  components  are  satisfied. 

A skeleton  template  of  the  Saber  LIGBT  model  is  shown  in 
Fig. 5.1  where  each  section  performs  the  following  functions:  The 
LIGBT  header  defines  the  A (anode),  G (gate),  and  K (cathode)  termi- 
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Template  LIGBT  A G K = model 
electrical  A,G,K  # node  type 
struct! 
number 

# user  input  parameters,  and  the  parameter  values 
be  changed  in  the  Saber  net-Hst. 

# local  declaration 

} 

{ 

parameters  { 

# parameters  caloilated  prior  to  simulation. 

} 

vEilues  { 

# values  calculated  as  a function  of  system  variables. 

} 

control  { 

# simulator  dependent  control  statements. 

} 

equations  { 

# equations  for  system  variables. 

} 


Fig.5.1  Skeleton  template  of  Saber  LIGBT  model. 
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nal  connection  points.  The  number  section  defines  the  user  input 
parameters,  and  the  parameters  can  be  changed  in  the  Saber  net-list. 
The  local  declarations  define  constants,  designate  internal  nodes,  and 
explicitly  define  the  additional  system  variables  needed  to  describe  the 
state  of  the  device  (usually  define  the  sample  point  of  the  variables  at 
each  interval).  The  parameters  section  is  used  to  calculate  quantities 
that  only  need  to  be  calculated  once  at  the  beginning  of  the  simulation. 
Quantities  that  are  functions  of  the  system  variables  (explained  in 
Chapter  3)  are  implemented  in  the  value  section.  The  control  section 
contains  information  about  nonlinear  model  relationships  and  com- 
mands to  aid  convergence.  Finally,  the  equations  section  describes  how 
the  quantities  calculated  in  the  values  section  are  assembled  to  solve 
for  system  variables. 

5.2.2  LIGBT  Model  Formulation 

To  implement  the  LIGBT  model  presented  in  Chapter  3 into 
the  Saber  circuit  simulator,  the  model  is  formulated  such  that  the  cur- 
rent between  each  of  the  terminal  nodes  are  expressed  in  terms  of  the 
nonlinear  functions  of  the  system  variables  and  in  terms  of  the  time 
rate  of  change  of  these  functions  of  the  system  variables.  The  expres- 
sions in  Chapter  3 are  implemented  in  the  values  section  of  the  Saber 
template,  and  these  values  are  used  in  the  equations  section  of  the 
template  to  describe  the  interconnection  of  the  components  of  current 
through  each  element  in  Fig.  5.1.  Figure  5.2  shows  a detailed  LIGBT 
equivalent  circuit  superimposed  on  a schematic  of  the  structure  of  one 
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(K)  (G)  (A) 


Cathode  Gate  Anode 


Fig.5.2  Network  representation  of  the  SABER  LIGBT  model 
superimposed  on  schematic  of  structure. 
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of  the  many  thousand  cells  of  an  n-channel  LIGBT.  The  elements  of  the 
circuit  of  Fig.  5.1  represent  the  nonlinear  physical  phenomena  associ- 
ated with  each  region  of  the  device  structure.  The  basic  LIGBT 
equivalent  circuit  of  a lateral  PNP  BJT/LIGBT  subcell  is  supplied  base 
current  by  the  LDMOST/LIGBT  subcell.  The  other  components  con- 
nected between  the  emitter  (E),  base  (B),  and  collector  (C)  node  are 
associated  with  the  lateral  PNP  BJT/LIGBT  subcell,  and  those  con- 
nected between  gate  (G),  source  (S),  and  drain  (D)  nodes  are  associated 
with  the  LDMOST/LIGBT  subcell. 

The  components  associated  with  steady-state  characteristics 
of  the  LIGBT  model  are  I^^,  Vepi,  Vcisc>  Imos»  potential 

drops  due  to  drift  and  diffusion  are  distributed  throughout  the  base 
region,  and  the  drift  terms  of  the  ambipolar  transport  equations  are 
coupled  [HEF88].  Thus,  both  the  base  and  collector  components  of  cur- 
rent contribute  to  the  resistive  potential  drop  Vj;pj.  V^jg^  represents 
the  potential  drop  across  current-induced-space-charge  region.  The 
current  components  of  I^,  Ij^,  and  Imos  represent  anode,  collector,  and 
LDMOST  currents,  respectively. 

The  components  associated  with  dynamic  characteristics  can 
be  divided  into  two  groups:  1)  the  lateral  PNP  BJT/LIGBT  subcell  com- 
ponents, and  2)  LDMOST/LIGBT  subcell  components.  The  components, 

^9^,  and  are  associated  with  the  lateral  PNP  BJT/LIGBT 
at  dt  dt 

subcell,  and  the  components,  — — j^,  and  — are  associated 

at  dt  dt 

with  LDMOST/LIGBT  subcell. 
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5.2.3  Saber  Implementation 

The  expressions  in  Chapter  3 are  implemented  in  the  values 
section  of  the  Saber  LIGBT  model  template  and  are  functions  of  the 
light  system  variables,  the  LIGBT  device  and  model  parameters  (see 
Chapter  4),  and  the  physical  constants  of  silicon.  The  system  variables 
for  the  LIGBT  are  the  node  voltages  and  the  explicitly  defined  system 
variables:  I^pi  and  Icisc*  The  first  seven  equations  in  the  values  sec- 
tion of  LIGBT  model  template  (see  Appendix)  evaluate  the  voltage 
differences.  The  notation  V(a)  is  the  MAST  syntax  for  the  voltage  at 
node  a.  The  quantities  evaluated  in  the  values  section  are  used  in  the 
equations  shown  in  Fig.  5.3  to  describe  the  currents  through  each  of 
the  elements  of  Fig. 5. 2,  and  to  describe  the  expressions  that  define  the 
explicitly  defined  system  variables:  I^pj  and  Icisc* 

The  first  twelve  statements  in  the  equations  section  of  Fig. 

5.3  describe  the  currents  between  node  pairs  of  Fig.5.2  in  terms  of  the 
system  variables.  The  notation  i(a->b)  indicates  that  a component  of 
current,  given  by  the  expression  on  the  right-hand  side  of  +=,  flows 
from  node  a to  node  b.  The  currents  through  the  gate-source  and  drain- 
source  capacitances  are  formulated  as  the  time  derivative  of  the 
charges  Qqs  and  Qds,  where  d_by_dt  is  the  Saber  time  derivative 
operator,  whereas  the  current  components  having  no  time  derivatives 
such  as  Imos>  Icisc»  Ia»  and  Ik  are  evaluated  in  the  values  section. 

The  last  two  statements  in  the  equations  section  give  the 
expressions  that  must  be  satisfied  to  determine  each  of  the  explicitly 
defined  system  variables.  The  notation  of  these  statements  indicates 
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equations  { 

i(G->B) 

+=d_by_dt(QGD) 

i(G->K) 

+=d_by_dt(QGS) 

i(S->K) 

+=IMOS 

i(B->S) 

+=ICISC 

i(B->K) 

+=d_by_dt(QCISC) 

i(A->E) 

+=IA 

i(E->B) 

+=IEPI 

i(A->B) 

+=d_by_dt(QBE) 

i(C->B) 

+=d_by_dt(QBC) 

i(B->C) 

+=d_by_dt(QJK) 

i(B->C) 

+=IK 

i(C->K) 

+=VCK/RP 

lEPI; 

v(E)-v(B)=VEPI 

ICISC: 

v(B)-v(S)=VCISC 

} 


Fig.5.3 


Equations  section  of  Saber  LIGBT  model  template. 
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that  the  system  variable  on  the  left-hand  side  of  the  colon  (:)  is  to  be 
solved  for  by  the  simulator  such  that  the  expression  on  the  right-hand 
side  of  the  colon  is  satisfied.  That  is,  in  addition  to  iterating  the  node 
voltages  until  Kirchhoff’s  current  law  is  satisfied  at  each  node.  The 
simulator  also  iterates  the  explicitly  defined  system  variables  until 
each  of  the  expressions  on  the  right-hand  sides  of  the  colons  are  satis- 
fied. System  variables  are  introduced  in  this  manner  to  describe 
components  of  current  that  cannot  be  expressed  as  explicit  functions  of 
the  node  voltages.  System  variables  are  also  introduced  to  describe  the 
time  derivatives  of  capacitor  voltages  where  the  capacitance  formula 
cannot  be  integrated  to  obtain  an  expression  for  the  capacitor  charge 
and  hence  the  nonlinearities  cannot  be  included  within  the  argument 
of  the  Saber  time  derivative  operator. 

In  the  LIGHT  template,  the  system  variables  Igpi  and  Icisc 
are  introduced  so  that  general  expressions  for  the  voltage  sources  can 
be  readily  implemented. 

5.2.4  Techniques  Used  to  Ensure  Convergence 

The  nonlinear  solution  algorithm  (iteration  algorithm)  of  the 
Saber  circuit  simulator  is  unique  in  that  all  of  the  nonlinear  expres- 
sions are  evaluated  in  conjunction  with  an  array  of  sample  points  for 
each  of  the  independent  variables.  At  the  sample  points,  the  nonlinear 
expressions  are  fully  evaluated,  while  at  intermediate  values  of  the 
independent  variables,  the  nonlinear  functions  are  evaluated  using 
multidimensional  linear  interpolations.  One  of  the  benefits  of  this 
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algorithm  in  implementing  model  equations  is  that  the  partial  deriva- 
tives of  the  model  equations  with  respect  to  the  system  variables  are 
not  required  in  order  for  the  simulator  to  iterate  the  system  variables 
in  a manner  that  converges  to  the  solution  of  the  system  of  nonlinear 
equations.  Thus,  the  model  equations  can  be  implemented  directly  in  a 
straightforward  manner  as  describe  above. 

The  range  and  the  density  of  the  sample  points  can  be  tai- 
lored to  the  nature  of  the  model  nonlinearities  by  specifying  the 
sample  points  arrays  in  the  control  section  of  the  Saber  template, 
which  is  shown  in  Fig.5.4.  The  simulator  also  provides  the  capability 
to  control  the  maximum  step  size  that  a given  variable  can  take 
between  successive  iterations  (Newton  step)  in  regions  where  the  non- 
linear model  functions  have  discontinuous  partial  derivatives  with 
respect  to  the  system  variables.  This  feature  of  the  Saber  simulator  is 
beneficial  in  implementing  models  which  use  different  expressions  to 
describe  different  regions  of  operation.  For  example,  Imos  is  described 
by  a different  expression  for  gate  voltages  above  and  below  Vgs=Vth. 
Newton  steps  for  the  independent  variables  can  be  introduced  near  the 
transition  between  different  regions  by  specifying  the  Newton  step 
arrays  in  the  control  section.  The  Newton  steps  tend  to  confine  the 
iterations  of  the  independent  variables  to  the  regions  where  they  are 
introduced  so  that  the  variables  do  not  overshoot  the  transition 
regions  during  iterations. 

In  implementing  the  LIGBT  model  into  the  Saber  circuit  sim- 
ulator, the  equations  in  Chapter  3 must  be  formulated  such  that  they 
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struc  sa_point  {number  bp,inc;}\ 

sva[*]=[(-100000, 1.0), (0,0.001), (0.5, 0.001), (1,0.1), (100000,0)], 

sveb[*]=[(-100000,0.01),(0,0.001),(0.5,0.001),(50000.0,0)], 

svch[*]=[(-50000,10),(0,0.1),(100000,0)j, 

svbk[*]=[(-50000,10),(0,.01),(1.0,1.0),(100000,0)], 

nv[*]=[(4.5,0.01),(5.5,0)l 


Control_section  { 

sample_points(VAK,sva) 

sample_points(VEB,sveb) 

sample_points(VBK,svbk) 

sample_points(VCH,svch) 

newton_step(VGS,nv) 


} 


Fig.5.4  Sample  points  and  control  section  of  Saber  LIGBT  model 
template. 
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are  continuous  and  nonsingular  in  the  range  that  the  system  variables 
may  take  during  iterations.  For  example,  the  expression  for  Xq  (base- 
collector  depletion  region  width)  is  not  valid  for  Vbk  < -0.6V,  and 
although  no  physical  solutions  exist  in  this  range,  it  is  necessary  to 
provide  a value  in  this  range,  because  the  base-collector  voltage  can 
enter  the  VBK  < -0.6V  range  during  the  iterations  necessary  to  find 
the  physical  solution.  However,  the  values  of  the  expressions  in  the 
range  where  there  are  no  solutions  are  only  important  during  itera- 
tions and  have  no  effect  on  the  results  of  the  model.  Therefore,  the 
functions  are  given  values  outside  the  physical  range  of  operation  such 
that  they  are  continuous  at  some  point  arbitrarily  near  the  singular- 
ity,  if  it  exists,  and  which  are  well  defined  beyond  the  range  of  physical 
operation. 

Model  Evaluation  and  Validation 

The  Saber  LIGBT  model  described  in  the  last  section  is  evalu- 
ated using  various  test  circuits  to  examine  the  device  behavior  for  the 
range  of  static  and  dynamic  conditions  in  which  the  device  is  intended 
to  be  operated.  We  have  four  general  external  circuits  to  test  the  Saber 
LIGBT  model:  1)  a series  resistor-inductor  load  and  resistive  gate 
drive  circuit,  2)  flyback  DC-to-DC  converter,  3)  soft-switching  circuit 
with  resonant  inductor  and  capacitor,  and  4)  full  bridge  inverter 
employing  four  LIGBT’s. 


101 


5.3.1  A Senes  Realstor-Inductor  Load  and  Gate  Drive  Circuit 

A series  resistor-inductor  load  and  resistive  gate  drive  test  cir- 
cuit, shown  in  Fig.  5.5,  is  used  to  verify  the  transient  model.  The  load 
circuit  state  equation  for  the  circuit  of  Fig.5.5  is 

= (5.1) 

where  /,.  = for  this  circuit.  The  gate  current  for  the  circuit  of  Fig. 
3.12  is  given  by 


gs' 


where  the  gate  pulse  generator  voltage  is  given  by 


0 

^gon  ~ ^^rise 


gon 

^gon  ^fall  ~ ^hall  ^nff+  ^fnll) 

0 


(5.2) 


(5.3) 


The  pulse  generator  rise  and  fall  times  are  described  by  the  second 
and  fourth  cases  in  Eqn.(5.3). 

The  corresponding  net-list  is  shown  in  Fig.  5.6.  The  state- 
ments of  Fig.5.6  describe  each  of  the  six  circuit  elements  of  Fig.  5.5, 
where  the  name  of  the  given  template  that  is  used  to  model  the  compo- 
nent is  given  on  the  left-hand  side  of  the  period  (.)  and  the  name  of 
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Fig.5.5 


Circuit  configuration  of  an  LIGBT  with  a series  resistor- 
inductor  load  and  a resistive  gate  drive. 
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LIGBT.ml  a g 0=  model=(dvw=0.1,tauh=2.6e-6) 

r.rl  al  a =30 

1.11  aa  al  =80u 
r.rg  ggg  =0.2K 

v.vg  gg  0 = tran  =(pulse=(0,15,0.5u,20n,20n,20u,40u)) 
v.va  aa  0 = dc  =300.0 


Fig. 5.6  Saber  net-list  for  a series  resistor-inductor  load  and  a 
resistive  gate  drive  circuit. 
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specific  component  within  the  circuit  is  given  on  the  right-hand  side  of 
the  period.  The  resistor,  inductor,  and  voltage  supply  templates  are 
provided  within  the  Saber  template  library,  and  the  user-defined 
template  is  described  in  Section  5.2.  The  remaining  columns  to 
the  left  of  the  equal  sign  define  the  connection  points  for  each  compo- 
nent, where  the  names  of  these  connection  points  are  defined  in  Fig. 
5.5.  The  parameters  of  each  component  that  are  to  be  changed  from 
the  default  values  of  the  generic  templates  are  given  to  the  right  of  the 
equal  signs. 

Figure  5.7  shows  the  anode  voltage  waveform  of  Saber  LIGHT 
model  for  different  values  of  gate  resistance.  The  value  of  the  voltage 
overshoot  at  turn-off  for  a stiff  gate  drive  (Rg  = 0.2  KG)  validates  the 
implementation  of  base-collector  depletion  charge  Qjc  which  domi- 
nates the  effective  output  capacitance.  The  turn-off  delay  time  for  the 
different  values  of  gate  resistances  validates  the  implementation  of 
the  low-voltage  gate-drain  capacitance.  The  anode  voltage  overshoot 
for  different  gate  resistances  also  validates  the  high-voltage  gate- 
drain  capacitance.  Figure  5.8  shows  the  anode  current  waveform  for 
three  different  values  of  gate  resistance.  The  slowly  deca)dng  portion 
of  the  turn-off  current  waveforms  validates  the  implementation  of 
anode  current  1^  and  non  quasi-static  charges  Qbe  and  Qgc.  The  gate 
voltage  and  gate  current  waveforms  for  different  values  of  gate  resis- 
tance is  shown  in  Fig.  5.9  and  Fig.  5.10,  respectively. 

These  characteristics  are  obtained  using  a Saber  command 
file  that  allows  you  to  perform  a parameter  sweep.  In  this  parameter 
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Fig. 5. 7 Saber  simulated  anode  voltage  waveforms  for  three  dif- 

ferent gate  resistances  (Rg  = 0.2K,  2.2K,  and  4.2K). 


Anode  Current  (amps) 
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Fig. 5. 8 Saber  simulated  anode  current  waveforms  for  three  dif- 

ferent gate  resistances  (Rg  = 0.2K,  2.2K,  and  4.2K). 
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Fig.5.9  Saber  simulated  gate  voltage  waveforms  for  three  dif- 
ferent gate  resistances  (Rg  = 0.2K,  2.2K,  and  4.2K). 
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Fig.5.10  Saber  simulated  gate  current  waveforms  for  three  dif- 
ferent gate  resistances  (Rg  = 0.2K,  2.2K,  and  4.2K). 
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sweep,  we  use  a range  of  gate  resistance  from  0.2K  ohms  to  4.2K  ohms, 
with  a step  size  of  2k  ohms.  This  means  that  there  is  one  execution  of 
the  transient  analyses  for  each  of  the  gate  resistance  values  0.2K, 
2.2K,  and  4.2K. 

5.3.2  Flyback  DC-to-DC  Converter 

The  circuit  of  flyback  DC-to-DC  converter  is  shown  in  Fig. 
5.11.  When  the  LIGBT  Q turns  on,  the  diode  D becomes  reverse-biased 
and  the  energy  flows  into  inductor  core  L.  When  Q turns  off,  the 
energy  stored  in  the  inductor  core  L causes  the  current  flow  through 
the  diode  in  the  counter-clockwise  direction.  In  this  flyback  converter, 
the  voltage-transfer  ratio  is 


^ouf  _ —D 

^in  ■ ra 


(5.4) 


where  D = to^/Tg  is  the  switch  duty  ratio.  Eqn.(5.4)  shows  that  the  volt- 
age transfer  ratio  in  a flyback  converter  depends  on  switch  duty  ratio 
D.  In  this  simulation,  D is  set  to  1/2  so  that  the  voltage  transfer  ratio 
becomes  -1. 

Figure  5.12  shows  the  anode-cathode  voltage  during  the 
operation  of  the  converter  circuit.  When  the  gate  voltage  is  on,  V^k  is 
almost  zero  (usually  less  than  2%  or  3%  of  the  peak  voltage),  and  when 
the  gate  voltage  is  off,  approaches  nearly  240  volts.  As  a switching 
device,  it  is  desirable  that  the  LIGBT  device  has  low  on-state  voltage 
and  high  off-state  voltage  across  the  anode  and  cathode  terminals. 
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ton  ~ 500n 


D = 


Tg  = ln 


^oi/Tg 


Q (LIGBT)  ^ (diode) 


Fig. 5. 11  The  flyback  DC  to  DC  converter  circuit.  The  duty  cycle  is 
set  1/2  to  make  output  voltage  of  -Vjn. 
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Fig. 5. 12  The  anode-cathode  voltage  during  on-off  condition  of  the 
gate  voltage. 
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Fig.5.13  Output  voltage  of  the  flyback  DC  to  DC  converter.  The 
ripple  can  be  eliminated  by  increasing  the  value  of 
capacitor. 
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The  output  voltage  across  the  resistor  R is  shown  in  Fig. 5. 13. 
Two  distinctive  behavior  is  superimposed  in  this  figure,  the  high  fre- 
quency behavior  by  a switching  circuit  and  low-frequency  behavior  of 
RLC  circuit  responding  to  step  input.  The  high-frequency  ripple  can  be 
reduced  by  increasing  the  value  of  capacitor  C.  The  output  voltage  has 
a initial  transient  condition  and  settles  down  to  steady-state  of  -100 
volts. 

5.3.3  Soft-Switching  Circuit  with  LC 

The  application  of  soft-switching  inverters  is  being  actively 
considered  by  many  manufacturers.  Most  available  IGBT’s  are  pres- 
ently designed  for  hard-switching  applications,  very  little  data  is 
available  in  the  literature  on  device  behavior  under  soft- switching 
conditions.  The  Saber  simulation  of  soft- switching  device  has  been 
motivated  by  a need  to  understand  IGBT  operation  under  the  different 
operating  conditions. 

Figure  5.14  shows  the  soft- switching  circuit  with  resonant 
inductor  and  capacitor.  The  test  soft-switching  circuit  used  for  Saber 
simulation  has  been  obtained  from  Kurnia  et  al.  [KUR92].  The  initial 
DC  condition  of  the  soft-switching  circuit  is  Vout  = Vin.  When  gate  volt- 
age is  applied  with  a step  function,  the  LIGBT  turns  on  and  the 
current  in  inductor  increases.  Therefore,  there  is  no  current  charging 
the  capacitor  C.  The  inductor  current  is  sinusoidal,  but  it  is  distorted 
when  the  LIGBT  switching  device  is  on.  When  the  LIGBT  turns  off,  the 
current  in  inductor  flows  to  charge  the  capacitor,  and  the  current 
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Fig.5.14  Soft-swtiching  circuit  with  resonant  inductor  and 
capacitor. 
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Fig. 5. 15  Saber  simulated  transient  capacitor  voltage,  gate  volt- 
age, and  inductor  current  waveforms. 
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waveform  becomes  sinusoidal.  When  the  voltage  across  the  capacitor 
becomes  higher  than  the  input  voltage,  the  inductor  current  flows  in 
the  reverse  way,  which  means  the  inductor  current  becomes  more  and 
more  negative.  Eventually,  the  capacitor  voltage  becomes  zero  as  the 
inductor  discharges  the  capacitor.  Even  when  the  voltage  across  the 
capacitor  becomes  zero,  the  inductor  current  still  continues  to  flow  its 
reverse  way.  This  makes  the  parasitic  diode  D turn  on,  and  the  diode 
on-state  voltage  drop  causes  the  negative  voltage  in  capacitor. 

When  the  inductor  current  approaches  zero,  the  gate  pulse 
generator  is  set  to  turn  the  gate  voltage  on  to  avoid  hard-switching, 
and  the  whole  cycle  repeats. 

5.3.4  Full-Bridge  Inverter 

A full-bridge  inverter  employing  four  LIGBTs  is  shown  in 
Fig. 5. 16.  This  inverter  consists  of  two  one-leg  inverters  and  is  pre- 
ferred over  other  arrangements  in  higher  power  ratings.  With  the 
same  DC  input  voltage,  the  maximum  output  voltage  of  the  full  bridge 
inverter  is  twice  that  of  the  half  bridge  inverter.  This  implies  that  for 
the  same  power,  the  output  current  and  the  switch  currents  are  one- 
half  of  those  for  a half  bridge  inverter. 

In  Fig.5.16,  the  diagonally  opposite  switches  (LIGBT.Ml, 
LIGBT.M4)  and  (LIGBT.M2,  LIGBT.M3)  from  the  two  legs  are  switched 
as  pairs  1 and  2,  respectively.  Figure  5.17  shows  the  turn-off  transient 
waveforms  of  the  current  and  voltage  for  the  LIGBT.M2.  The  current 
remains  steady  while  the  voltage  across  the  LIGBT.M2  switching 
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Fig.5.16  Schematic  circuit  of  full  bridge  inverter  employing  four 
LIGBTs. 
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device  rises.  Only  after  the  voltage  approaches  its  off-state  value  does 
current  begin  to  fall.  These  are  standard  inductive  switching  wave- 
forms as  commonly  seen  in  switch-mode  power  converters. 


5.4  Summary 

The  physics-based  model  for  the  LIGHT  has  been  imple- 
mented into  the  Saber  circuit  simulator.  The  techniques  for 
implementing  the  LIGHT  model  equations  into  the  Saber  circuit  simu- 
lator are  essential  for  accurate  simulation  of  other  power 
semiconductor  devices.  The  model  has  been  verified  for  various  exter- 
nal circuit  conditions  and  for  the  full  range  of  static  and  dynamic 
conditions  in  which  the  LIGHT  is  intended  to  be  operated.  The  model 
has  also  been  shown  to  be  suitable  for  simulating  the  behavior  of  the 
LIGHT  for  general  purpose  external  circuit  conditions  and  designing 
the  various  power  circuits. 


CHAPTER  6 

CONCLUSION  AND  SUGGESTIONS  FOR  FUTURE  WORK 


In  this  dissertation,  an  accurate  charge-based  model  for  the 
LIGHT  has  been  developed  and  shown  to  be  useful  for  optimal  device/ 
circuit  CAD  of  IGBTs.  Improvements  on  both  subcell  models  of  the 
LIGHT  have  resulted  in  accurate  simulation  results  without  parame- 
ter optimization.  In  the  LDMOST/LIGBT  subcell,  we  have  modeled  the 
current-induced-space-charge  (CISC)  region.  The  model  of  CISC  region 
has  been  shown  to  be  of  great  importance  for  improving  the  accuracy  of 
the  LDMOST  model.  In  the  lateral  PNP  BJT/LIGBT  subcell,  we  have 
derived  analytical  geometrical  factors  based  on  the  ambipolar  trans- 
port equation.  The  lateral  and  quasi-vertical  geometrical  factors  are 
needed  in  the  two-dimensional  anode  current  and  thus  in  accurate 
modeling  for  steady-state  and  transient  conditions  for  the  LIGHT.  The 
inclusion  of  the  Kirk  effect  in  our  LIGHT  model  also  improved  the 
accuracy  of  steady-state  characteristics  especially  for  LIGHT  with 
short  distance  between  the  anode  and  the  cathode.  The  addition  of  Qbc 
term  in  the  LIGHT  transient  model  also  improved  the  accuracy  of  tran- 
sient current  and  voltage  waveforms  of  the  LIGHT  for  general  loading 
conditions. 

The  physical  insights  gained  from  two-dimensional  numerical 
simulations  and  experimental  measurements  of  test  devices  were  used 
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to  develop  our  charge-based  LIGBT  models.  The  generalized  methodol- 
ogy for  physics-based  modeling  of  IGBTs  based  on  the  two-dimensional 
numerical  device  simulator  PISCES  can  be  useful  to  model  the  power 
devices  with  wide  and  lightly  doped  regions.The  LIGBT  model  was 
implemented  in  the  Saber  circuit  simulator,  by  which  the  model  equa- 
tions are  solved  semi-numerically  within  the  nodal  analysis 
framework  of  Saber.  Saber  is  a very  powerful  modeling  tool  for  power 
devices  which  have  widely  varying  cell  structures  even  within  one  cat- 
egory of  device. 

Our  new  model  was  verified  by  measurements  of  specially 
designed  test  devices,  using  simulations  with  model  parameters  evalu- 
ated from  device  structural  information  and  static  and  dynamic 
measurements. 

Parameter  extraction  method  is  developed.  We,  however,  used 
PISCES  simulation  results  due  to  unavailability  of  accurate  measure- 
ment data.  The  model  parameters  are  obtained  sequentially  by 
selecting  features  of  the  device  characteristic  that  isolate  parameters, 
and  using  the  parameters  obtained  from  previous  extraction  steps  to 
calculate  the  model  parameters  from  the  next  PISCES  simulated 
characteristics. 

The  Saber  LIGBT  model  has  been  used  to  simulate  various 
external  circuits  and  has  been  proven  valid  for  simulating  and  design- 
ing the  power  electronic  circuits. 

Based  on  the  research  discussed  herein,  we  suggest  the  fol- 
lowing topics  for  future  research  consideration. 
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First,  we  recommend  a study  of  regional  models  for  IGBT. 
Each  regional  building  block  represents  the  physics-based  model  of 
each  basic  semiconductor  region  in  the  IGBT.  Then  IGBT  model 
becomes  a simple  association  of  several  of  the  basic  semiconductor 
region  models.  These  regional  building  blocks  can  be  reused  in  model- 
ing of  other  power  devices  with  different  structures  such  as  PIN  diode 
and  power  transistor. 

Second,  we  recommend  the  refinements  of  LDMOST  current 
model.  The  first-level  of  MOSFET  current  equation  has  been  used  in 
our  LIGBT/LDMOST  subcell  model.  More  accurate  description  is 
needed  in  the  model  of  the  LDMOST/LIGBT  subcell,  for  example  to 
account  for  secondary  and  parasitic  effects  in  the  LDMOST  channel 
region. 

Third,  we  recommend  the  implementation  of  our  LIGBT 
model  directly  into  SPICE  source  code.  The  direct  SPICE  implementa- 
tion requires  the  model  routine  to  evaluate  all  the  element  values  from 
the  device  terminal  voltages. 

Fourth,  we  recommend  a study  of  the  systematic  parameter 
extraction  with  empiricism.  The  parameter  extraction  can  be  auto- 
mated by  linking  ICCAP  to  Saber  circuit  simulator. 

Fifth,  we  recommend  the  comparison  between  the  measure- 
ment data  of  external  test  circuits  shown  in  Chapter  5 and  the  Saber 
simulation  results.  This  can  be  done  either  by  obtaining  actual  multi- 
cell LIGBT  devices  or  by  modifying  element  values  of  the  test  circuits 
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based  on  the  maximum  current  capacity  of  the  single-cell  LIGHT  V 
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APPENDIX  A 

LIGBT  PISCES  INPUT  FILE 


In  this  appendix,  the  PISCES  input  file  for  the  LIGBT  device 
provided  by  AT&T  Bell  Laboratories  is  listed.  The  PISCES  input  file 
details  the  geometry  input  and  processing  information  of  the  AT&T 
LIGBT  device. 
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PISCES  INPUT  FILE 

TITLE  LIGHT  PISCES  Simulation  Input  File 
+ : Simulation  Structure  Generation 


$ Mesh  Generation 

MESH  rectangu  jt.min=Q.Q  x.max=55.Q  v.min=Q.Q  v.max=2Q.Q 

+ smooth. k=l 

X.MESH  x.left=0.0  x.right=8.5  hl=0.5  ratio=1.0 

X.MESH  x.left=8.5  x.right=12.0  hl=0.25  ratio=1.0 

X.MESH  x.left=11.0  x.right=39.25  hl=2.0  ratio=1.0 

X.MESH  x.left=39.25  x.right=42.75  hl=0.25  ratio=1.0 

X. MESH  x.left=42.75  x.right=55.0  hl=1.0  ratio=1.0 

Y. MESH  node=l  location=-0.1 

Y.MESH  node=4  location=0.0 

Y.MESH  y.top=0.0  y.bottom=6.0  hl=0.25  ratio=1.0 

Y.MESH  y.top=6.0  y.bottom=20.0  hl=2.0  ratio=1.0 


$ Mesh  Elimination 
ELIMINATE  columns  x.min=8.75 
+ y.max=20.0 

ELIMINATE  columns  x.min=8.75 
+ y.max=20.0 


x.max=11.75  y.min=6.0 
x.max=11.75  y.min=6.0 


ELIMINATE 

ELIMINATE 

+ 

ELIMINATE 

+ 


columns  x.min=0.0  x.max=8.5  y.min=6.0  y.max=20.0 
columns  x.min=39.25  x.max=42.75  y.min=6.0 
y.max=20.0 

columns  x.min=39.25  x.max=42.75  y.min=6.0 
y.max=20.0 


$ Region  Definition 

REGION  number=l  oxide  iy.min=l  iy.max=4 

REGION  number=2  silicon  x.min=0.0  x.max=55.0  iy.min=4 

+ y.max=20.0 


$ Electrodes  : #l=Anode  #2=gate  #3=Cathode 
ELECTRODE  number=l  x.min=45.0  x.max=55.0  iy.min=l 
+ iy.max=4 

ELECTRODE  number=2  x.min=8.0  x.max=14.5  top 
ELECTRODE  number=3  x.min=0.0  x.max=7.0  iy.min=l  iy.max=4 

$ Substrate  & Epi-layer  Doping 

PROFILE  n-type  n.peak=2.0el4  x.peak=Q.Q  uniform  y.peak=0.0 
+ depth=20.0  outfile=hdop3.2 


$ P'  Body  Doping 

PROFILE  p-type  x.peak=0.0  x.right=8.0  n.peak=6.Qel7 
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+ y.peak=0.4  y.juncti=4.0  xy.ratio=0.9 

$ Deep  Doping 

PROFILE  p-type  x.peak=0.0  x.right=3.5  n.peak=1.0el9 
+ y.peak=0.0  y.juncti=6.0  xy.ratio=0.8 

$ Source  Doping 

PROFILE  n-type  x.peak=5.05  x.right=8.2  n.peak=5.0el9 
+ y.peak=0.0  y.char=0.55  xy.ratio=0.65 

$ P^  Anode  Doping 

PROFILE  p-type  x.peak=45.0  x.right=55.0  n.peak=1.0el9 
+ y.peak=0.0  y.junction=6.0  xy.ratio=0.8 

$ Specification  of  Fixed  Charge 
INTERFACE  qf=5.0el0 


$ 2-D  Plot  of  Initial  Grid 

PLOT.2D  bounjunc  l.bound=2  l.elect=l  title="LIGBT  Init" 

+ x.min=0.0  x.max=16.0  y.min=-0.1  y.max=8.0 

+ scale  device=xterm  pause 


$ Regrid  on  Doping 

REGRID  doping  logarith  smooth. k=l  ignore=l  cos.angl=2.0 
+ ratio=2.0  dopfile=bdop3.2 

PLOT.2D  grid  title="LIGBT  Doping  Regrid" 


+ x.min=0.0  x.max=12.0  y,min=-0.1  y.max=8.0 

+ device=xterm  pause 


$ Gate  Contact  Parameter 
CONTACT  number=2  n.poly 


$ Physical  Model 

MODELS  conmob  fldmob  srfmob2 


$ Symbolic  Factorization,  Solve  & Regrid  on  Potential 
SYMBOLIC  carriers=0  gummel  min.degr  strip 
METHOD  iccg  damped 

SOLVE  vl=0.0  v2=0.0  v3=0.0 

REGRID  potentia  smooth. k=l  ignore=l  cos.angl=2.0 

+ ratio=0.2  max.leve=l  dopfile=hdop3.2  outfile=pot3.mesh2 

PLOT.2D  bounjunc  l.bound=2  l.elect=l  title="LIGBT  Junction 
+ Structure" 

+ scale  device=xterm 


STOP 


APPENDIX  B 

SABER  LIGBT  MODEL  TEMPLATE 


In  this  appendix,  we  overview  the  Saber  implementation  of  the 
composite  LIGBT  model  developed  in  chapter  3.  All  of  the  elements  in 
the  network  representation  of  the  LIGBT  model  shown  in  Fig.5.2  are 
implemented  in  Saber  circuit  simulator. 

The  source  code  of  the  LIGBT  Saber  template,  which  details  the 
model  routine  for  the  LIGBT  model  shown  in  Fig.5.2,  is  listed  on  the 
following  pages. 
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########################################################### 

# 

#THIS  WORK  WAS  SUPPORTED  BY  THE  SRC  UNDER  CONTRACT 

# 92-SJ-285.  COPYRIGHT,  1993,  UNIVERSITY  OF  FLORIDA. 

# ALL  RIGHTS  RESERVED. 

# FOR  FURTHER  INFORMATION,  CONTACT  ; 

# TAEHOON  KIM  AND  D.E.  BURK 

# E-MAIL  ADDRESS  : THK@TCAD.EE.UFL.EDU 

# OR  DEB@TCAD.EE.UFL.EDU 

# FEBRUARY,  1ST,  1995 

# 

########################################################### 


# ligbt.sin  LIGBT.sin  (Physical  LIGBT) 

# AH  parmaters  are  in  M.K.S.  Unit 

# Temperature  dependent  parameters  are  installed 


element  template  light  A G K = model 

#Header  declaration 

electrical  A,G,K 
struc  { 

################################################### 
# User  input  paramter 

################################################### 

number  T = 300.0 
number  PD  = 5.0e22 
number  TAUH  = 2.6e-6 
number  VTH  = 4.8 
number  NEPI  = 6.0e20 
number  AD  = 1.0e25 
number  XJ  = 6.4e-6 
number  WCH  = 2.0e-6 
number  DVH  = 16.0e-6 
number  PA  = lO.Oe-6 
number  ALJ  = 8.0e-6 
number  WBV  = 65.0e-6 
number  WB  = 50.0e-6 
number  TOX  = l.Oe-7 
number  DVW  = 120.0e-6 
number  LSOL  = l.Oe-6 


129 


number  LDOL  = 6.0e-6 
number  KPLIN  = 0.96e-4 
number  KPSAT  = 0.41e-4 
} model  = 0 

{ 

######################################################## 
# Physical  constants 

######################################################## 

number  QUE  = 1.6022e-19  # Element  charge 
number  PI  = 3.14159 

number  BOLT=  1.38e-23  # Boltzman  constant 

number  EPSO=  3.45e-l  # Oxide  dielectric  constant 

number  EPSS=  l.Oe-10  # Silicon  dielectric  constant 

number  VS  = 1.07e5  # Saturation  velocity 


#Values  and  variables  declaration 


number  XKT,XP,REFT,PHIB,DNOMX,DNOMY,BR,COX,DN,DP, 
RTO,EG,NI,UN,UP,LA,TAUN,DA,JNO,RP,ARG,XLN, 
AREA, UNX,UPX,NEPIX,VT,WEPI, ADD, DNOMZ,UNZ, 
UNA,DNA,J1,MUC,TERM 

val  V VAK,VEB,VBK,VCK,VEPIX,VEPIY,VEPIZ,VC,VEPI, 
VCISC,VCH,VGS,VGD,VBKL,VSBE 
val  i IA,IMOS,IB,IK,IAL,IAO 

val  q QGS,QGD,QBL,QBV,QJK,QCISC,QLAU,QLAD,QLA, 

QLBUL,QLBUR,QLBD,QLB,QBCL,QBEL,QBCV, 
QBEV,QBE,QBC 
val  meter  XD,W,WX,XWM,WM 

val  nu  PO,XH,JAA,JAB,JA,JBX,JB,JMOS,JC,NCISC,XFAC, 

COFA,COFBA,QXT,AMOS,COFB,COFCA,COFCB,COFCBA, 
COFCD,COFCC,COFC,XINTA,XINTB,XINT,FGAL,CA,CU, 
CENY,MAXR,MINR,USEC,DSEC,PFGAV,SFGAV,FGAV, 
ACISC,CQL,PIC,  QXTT,QBCLE,JK,XD1,BETA,LLA, 
ALPHA,JN,DPP,KPSAT,KPLIN 
var  i iepi,icisc 


group 

group 

group 


{VCISC,VEPI}  V 

{IA,IMOS,IB,IK}  I 

{QGS,QGD,QBC,QBE,QJK,QCISC}  Q 
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####################################################### 
# internal  nodes 

####################################################### 

electrical  E,B,S,C 


# 


struc  sa_point  {number  bp,inc;}\ 


sva[*]=[(-10000, 1.0), (0,0.01), (0.5, 0.01), (1,0.1), (10000,0)], 
sveb[*]=[(-10000,0.01),(0,0.001),(0.5,0.001),(5000.0,0)], 
svch[*]=[(-50000,10),(0,0.1),(100000,0)], 
svbk[*]=[(-50000,10),(0,.01),(1.0,1.0),(100000,0)], 
nv[*]=[(4.5,0.01), (5.5,0)] 

# 


PARAMETERS  { 

EG  = 1.16-(7.02e-4*model->T**2.)/(model->T+1108.0) 

VT  = model->T*8.61e-5 
XKT  = BOLTZ*model->T 

ARG  = -EG/(XKT+XKT)+(  1.115 1278409 1/(BOLTZ*600.)) 

XP  = LIMEXP(QUE*ARG) 

REFT  = model->T/300. 

NI  = 1.33el6*REFT**1.5*XP 

PHIB  = VT*LN(model->NEPI*5.0e21/NI**2.) 

NEPIX  = model->NEPl*1.0E-6 
ADD  = model->AD*1.0E-6 

DNOMX  = 1.0+(NEPIX/(1.26el7*REFT**2.4))*0.88*REFT**-0.146 

DNOMZ  = 1.0+(ADD/(1.26el7*REFT**2.4))*0.88*REFT**-0.146 

DNOMY  = 1.0+(NEPIX/(2.35el7*REFT**2.4))*0.88*REFT**-0.146 

UNX  = 88.0*REFT**-0.57+(7.4e8*model->T**-2.33)/DNOMX 

UNZ  = 88.0*REFT**-0.57+(7.4e8*model->T**-2.33)/DNOMZ 

UPX  = 54.3*REFT**-0.57+(1.36e8*model->T**-2.23)/DNOMY 

UN  = UNX*1.0E-4 

UNA  = UNZ*1.0E-4 

UP  = UPX*1.0E-4 

BR  = UN/UP 

COX  = EPSO/model->TOX 

DN  = XKT/QUE*UN 

DNA  =XKT/QUE*UNA 

DP  = XKT/QUE*UP 

DA  = 2*DN*DP/(DN+DP) 

LA  = SQRT(DA*model->TAUH) 
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RTO  = model->DVW/model->WCH 
MUC  = UN/(l.+1.0E-6*(VGS-model->VTH)) 

TERM  =MUC*RTO*COX 

TAUN  = l./((l.+1.0e-17*ADD)/3.5e-5+8.3e-32*ADD**2.) 

XLN  = (DNA*TAUN)**0.5 

JNO  = QUE*DNA**0.5/(TAUN**0.5*model->AD) 

AREA  = model->DVH*model->DVW 

WEPI  = model->DVH-model->XJ 

Rp  = l./(model->PD*UP*QUE*model->DVW) 

J1  = QUE*model->NEPI*VS 

} 


# 


values  { 

########################################################### 

# Node  voltages 

########################################################### 

VAK  =v(A)-v(K) 

VGS  =v(G)-v(K) 

VEB  =v(A)-v(E) 

VBK  =v(B)-v(K) 

VCH  =v(S)-v(K) 

VGD  =v(G)-v(B) 

VCK  =v(C)-v(K) 

########################################################### 

# Calculating  Anode  current 

########################################################### 

PO  =(NI**2./model->NEPI)*LIMEXP(VEBA^T) 

IF(VBK<=0.0)  { 

XD=0 

} 

ELSE! 

XDl=(2*EPSS*(PHIB+VBK)/(QUE*model->NEPI))**0.5 
IF(JK<=0.0)  { 

XD=XD1 

} 

ELSE! 

XD=XDl*(l.+JK/Jl)**-0.5 

} 

} 

W=model->WB-XD 
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JAA=JN0*(P0/NI)**2. 
JAB=(QUE*DP*P0/LA)*1.0/TANH(W/LA) 
JA=(BR+1)/BR*(JAB+JAA) 
ALPHA=1.-(2.*DN/(DN+DP)*(1.-1./C0SH(W/LA))) 
BETA=ALPHA/(  1- ALPHA) 

DPP=DP 


########################################################### 
# Calculating  VEPI 

########################################################### 

IF(IEPI<=0)  { 

VEPI=0 

} 

ELSE  { 

VSBE=VT*LN(model->NEPl**2./NI**2.) 

IF(VEB<=VSBE)  { 

VEPI=0.0 

} 

ELSE  { 

XH=model->NEPI*SINH(W/LA)/PO 

XWM=LA*LN(XH+(XH**2.+1.0)**0.5) 

VEPIY=(BR-l)/(BR+l)*VT*LN(PO/model->NEPI) 

WM=W-XWM 
IF(WM>=W)  { 

VEPIZ=0.0 

} 

ELSE  { 

VEPIX=JA*SINH(W/LA)*LA/(QUE*(UN+UP)*PO) 

VEPIX=JA*SINH(W/LA)*LA/(QUE*(UN+UP)*PO) 

VEPIY=(BR-l)/(BR+l)*VT*LN(PO/model->NEPI) 

VEPIW=(LN((TANH(W/(2.*LA)))) 

VEPIV=VEPIW-LN((TANH((W-WM)/(2.*LA))))) 

VEPIZ=VEPIX*VEPIW 

} 

VEPI=VEPIZ+VEPIY 

} 

} 


########################################################### 
# Calculating  JB 

########################################################### 

JBXX=(QUE*PO*LA/(model->TAUH*SINH(W/LA))) 


JBX=JBXX*(COSH(W/LA)-1.0) 

JB=JBX+JAA 

JK=JA-JB 


133 


########################################################### 
# Calculating  IMOS 

########################################################### 

VC=VGS-model->VTH 
IF(VC<=0.0)  { 

IMOS=0.0 

} 

ELSE  { 

IF(VCH<VC)  { 

KPLINN=model->KPLIN*(model->DVW/120e-6) 

KPSATT=model->KPSAT*(model->DVW/120e-6) 

IMOS=KPLINN*(VGS-model->VTH-VCH/2.)*VCH 

} 

ELSE  { 

IMOS=KPSATT*/2.*(VGS-model->VTH)**2. 

} 

} 


########################################################### 
# Calculating  VCISC 

########################################################### 


IF(ICISC<=0)  { 

VCISC=0 

[ 

ELSE  { 

JC=QUE*VS*model->NEPI 
IF(XD<=model->WCH)  { 

WX=0.0 

VCISC=0.0 

} 

ELSE  { 

WX=XD-model->WCH 

AMOS=0.5*PI*WX*model->DVW 

JMOS=IMOS/AMOS 

NCISC=JMOS/(QUE*VS) 

VCISC=0.5*QUE/EPSS*NCISC*WX**2.0 

} 
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} 


########################################################### 

# Two  dimensional  analysis  • Lateral  side 
########################################################### 

XFAC=LIMEXP(2*model->WB/LA) 

COFA=(XFAC+1.0)/(XFAC-1.0) 

COFBA=(XFAC+1.0)*XFAC*2.0*model->XJ/((XFAC-1.0)*LA) 

COFB=(2.0*model->XJ*XFAC/LA-COFBA)/(XFAC-1.0) 

COFCA=1.0/12.0*XFAC*2.0*model->XJ*(12.0*model->XJ-LA)/LA**2 

COFCBA=(12.0*model->XJ-LA)/((XFAC-1.0)*LA**2) 

COFCB=1.0/12.0*(XFAC+1.0)*XFAC*2.0*model->XJ*COFCBA 

COFCD=(2.0*model->WB/LA)**2 

COFCC=LIMEXP(COFCD)*4.0*model->XJ**2/((XFAC-1.0)**2*LA**2) 

COFC=(COFCA-COFCB+COFCC)/(XFAC-1.0) 

PIC=PI/2.0 

XINTA=COFA+COFB*PIC**2+COFC*PIC**4+24.0*COFC*(1.0+PIC) 

X1NTB=2.0*COFB*(1.0+PIC)+4.0*COFC*PIC**3+12.0*COFC*P1C**2 

X1NT=XINTA-XINTB 

FGALL=model->DVW*XINT/(model->DVW+model->PA) 
FGAL=FGALL*  l./TANH(model->WB/LA) 

########################################################### 

# Two  dimensional  analysis  - Vertical  side 
########################################################### 


CA=model->WB/model->XJ 

CU=WEPI/model->XJ 

CENY=((model->WBV/2.0)**2-WEPI**2)/(2.0*WEPl) 

MAXR=SQRT(CENY**2+(model->WBV/2.)**2.) 

M1NR=MAXR-WEPI 

USEC=ACOS(model->XJ/MAXR) 

DSEC=ACOS(model->XJ/MINR) 

PFGAV=CA/2.0*(USEC+USEC**3/9.0+2.0*USEC**5/75.0) 

SFGAV=CA/2.0*(DSEC+DSEC**3/9.0+2.0*DSEC**5/75.0) 

FGAV=ABS(PFGAV-SFGAV) 

########################################################### 
# Charging  Terms 

########################################################### 

QGD=(model->LDOL*model->DVW*COX)*VGD 

QGS=(model->LSOL*model->DVW*COX)*VGS 

ACISC=0.25*PI*WX**2.0 
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QCISC=0.0*QUE*ACISC*model->DVW*NCISC 

CQL=LIMEXP(W/LA) 

QLAU=0.5*CQL+1./(2.*CQL)-1 
QLAD=0.5*CQL- l./(2.*CQL) 

QLA=QLAU/QLAD 

QLBUL=(0.5*(CQL-model->XJ)/LA)-(0.5*model->XJ/(LA*CQL)) 

QLBR=(CQL-l.)*(0.5*CQL*model->XJ/LA+0.5*model->XJ/(LA*CQL)) 

QLBUR=QLBR/(CQL+1.) 

QLBD=0.5*CQL-(1./(2.*CQD) 

QLB=(QLBUL-QLBUR)/QLBD 

QXT=COSH(W/LA)-1.0 

QXTT=0.5*(CQL**2+1.)/(CQL*SINH(W/LA)M1./SINH(W/LA)) 

QBL=QUE*model->DVW*LA*PO*model->XJ*QXTT 

QBCLE=0.5*(CQL**2-1.)*LA/(CQL*W*SINH(W/LA))-(1./SINH(W/LA)) 

QBCL=QUE*model->DVW*LA*PO*model->XJ*QBCLE 

QBEL=QBL-QBCL 

QBV=QBL*(FGAV/FGAL) 

QBEV=QBV*(QBEL/QBL) 

QBCV=QBV-QBEV 

QBE=QBEL+QBEV 

QBC=QBCL+QBCV 

QJK=QUE*XD*AREA*model->NEPI 

########################################################### 
# Calculating  Anode  current 

########################################################### 

IAO=JA*(model->DVW+model->PA)*model->XJ 

IA=IAO*(FGAL+FGAV) 

IK=IA-JB*AREA 

JK=JA-JB 

IB=JB*AREA 


} 

# 

Control_section  { 

sample_points(VAK,sva) 

sample_points(VEB,sveb) 

sample_points(VBK,svbk) 

sample_points(VCH,svch) 

newton_step(VGS,nv) 


} 
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equations  { 

i(G->B)  +=d_by_dt(QGD) 
i(G->K)  +=d_by_dt(QGS) 
i(S->K)  +=IMOS 
i(B->S)  +=ICISC 
i(B->K)  +=d_by_dt(QCISC) 
i(A->E)  +=IA 
i(E->B)  +=IEPI 
i(A->B)  +=d_by_dt(QBE) 
i(C->B)  +=d_by_dt(QBC) 
i(B->C)  +=d_by_dt(QJK) 
i(B->C)  +=IK 
i(C->K)  +=VCK/RP 

lEPI:  v(E)-v(B)=VEPI 

ICISC:  v(B)-v(S)=VCISC 

} 

} 
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